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The  prospect  of  sustainable  management  of  natural  tropical  forests  for  their  timber 

and  non-timber  benefits  faces  important  biological  and  technical  challenges  of  defining 

the  conditions  necessary  for  sustainability,  as  well  as  economic  impediments  concerning 

opportunity  costs  and  appropriate  incentive  mechanisms.  In  this  study,  an  optimization 

model  of  forest  management  is  employed  to  examine  the  conditions  necessary  for 

sustainable  management  of  a  Bolivian  tropical  dry  forest.  Integration  of  a  forest  growth 

model  in  the  optimization  enables  an  examination  of  harvesting  impacts  on  forest 

structure  and  composition  in  addition  to  the  economic  returns  to  management.  Results  of 

the  optimization  analyses  are  then  employed  to  examine  the  efficiency  of  alternative 

fiscal  regulatory  mechanisms  in  promoting  more  sustainable  logging  behavior. 

Simulation  results  indicate  that  sustainable  management  for  timber  production  is  both 

feasible  and  profitable  in  the  dry  forest  site.  Alternative  regulatory  policies  may  achieve 

sustainability  goals  more  efficiently  than  current  Bolivian  forest  policy,  but  compromises 

between  the  efficiency  and  practicality  of  regulatory  approaches  are  warranted. 


CHAPTER  1 
INTRODUCTION 

The  success  of  initiatives  promoting  sustainable  forest  management  (SFM)  in 
natural  tropical  forests  will  depend  on  demonstration  that  it  is  not  only  technically 
feasible,  but  also  equally  decisive,  that  it  is  economically  viable.  Critical  to  the  effective 
conservation  of  tropical  forest  ecosystems  managed  for  timber  will  be  tangible  incentives 
for  improved,  sustainable  management  by  forest  landowners  and  the  logging  industry. 
The  management  of  natural  tropical  forests  for  timber  is  a  complex  and  contentious  issue. 
Its  complexity  derives  from  the  nature  of  the  mixed-species,  uneven-aged  tropical  forests 
themselves,  which  confound  attempts  to  estimate  the  conditions  necessary  for 
sustainable,  productive  management.  Its  contentiousness  derives  from  the  extraordinary 
nonmarket  value  of  the  biodiversity  and  ecological  service  flows  conserved  by  tropical 
forests,  the  increasing  scarcity  of  tropical  forests,  and  the  deleterious  consequences  of 
their  mismanagement. 

Although  the  technical  challenges  to  defining  and  implementing  sustainable  forest 
management  (SFM)  are  imposing,  most  authors  agree  that  the  principal  impediment  to  its 
greater  adoption  is  SFM's  financial  inferiority  relative  to  conventional  logging  practices, 
and  alternative  land  uses  (Verissimo  et  at.  1992,  Kishor  and  Constantino  1993,  Rice  et  al. 
1997,  Pearce  et  al.  2003).  Given  this  conclusion,  and  the  desirability  of  a  sustainable 
solution  to  unabated  demands  for  the  harvest  of  tropical  timbers,  there  is  urgent  need  for 
empirically  rigorous  examination  of  the  marginal  costs  of  more  sustainable  management, 
and  the  trade-offs  that  may  be  optimal  to  achieve  both  production  and  conservation  goals. 
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Recent  studies  of  Dipterocarp  forests  in  peninsular  Malaysia  examined  such  trade-offs 
between  the  financial  and  ecological  benefits  of  natural  forest  management  (Boscolo  and 
Buongiorno  1997,  Ingram  and  Buongiorno  1997,  and  Boscolo  and  Vincent  2000).  The 
present  research  contributes  to  earlier  work  by  examining  the  conditions  necessary  for, 
and  the  biological  consequences  and  economic  costs  of  sustainability,  of  a  forest  in  the 
eastern  Bolivian  lowlands. 

The  Bolivian  forestry  law  was  revised  in  1996  (Ley  1 700)  to  more  effectively 
promote  sustainable  management  of  the  nation's  forest  resources,  which  are  allocated  in 
renewable  concessions  to  the  logging  industry  (Art.  #1).  The  present  study  offers  a 
contribution  to  the  analysis  of  whether  or  not  new  forestry  legislation  is  likely  to  achieve 
these  sustainability  goals.  The  study  offers  an  economic  comparison  of  management 
alternatives  under  current  and  alternative  regulatory  constraints,  and  a  robust  analysis  of 
the  biological  impact  of  such  alternatives  for  stands  representative  of  the  Chiquitano 
tropical  dry  forests  of  eastern  Bolivia. 

The  study  focuses  on  800  hectares  (ha)  of  seasonally  dry  forest  located  near  the  Las 
Trancas  community  (16°13'S,  61°50'W)  in  the  Lomerio  region  south  of  Conception, 
Santa  Cruz,  Bolivia  The  forest  is  representative  of  the  Chiquitano  tropical  dry  forest 
ecosystem,  which  lies  in  a  transition  zone  between  the  humid  forests  of  the  Amazon 
Basin  and  the  thorn  scrub  of  the  Gran  Chaco  (Killeen  et  al.  1998).  Numerous  timber 
concessions  have  been  awarded  in  the  Chiquitano  forests,  which  are  valued  for  their 
stocks  of  tropical  hardwood  species  such  as  Spanish  cedar  (Cedrelafissilis),  Spanish  oak 
{Amburana  cearensis),  and  "morado"  (Machaerium  scleroxylon). 


Objectives  of  the  study  are  to  portray  the  optimal  behavior  of  a  forest 
concessionaire,  whose  goal  is  to  maximize  the  financial  returns  from  logging  and  whose 
management  actions  are  circumscribed  by  biological,  market,  and  regulatory  constraints. 
The  study  aims,  moreover,  to  determine  the  conditions  necessary  to  maximize  the  returns 
from  logging  while  meeting  the  goals  of  sustained  timber  yield  and  maintenance  of 
desired  forest  conditions.  The  study  comprises  three  principal  elements,  each  of  which  is 
fully  developed  in  the  ensuing  chapters.  These  elements  include: 

•  Estimation  of  a  forest  growth  and  yield  model  as  a  basis  for  establishing  the 
biological  conditions  by  which  management  is  constrained 

•  An  economic  optimization  comparing  the  profitability  and  forest  impacts  of 
alternative  management  regimes 

•  An  analysis  of  regulatory  policies  as  means  of  providing  incentives  to 
concessionaires  for  more  sustainable  forest  management  behavior 

Together  these  components  respond  to  fundamental  issues  of  sustainable  forest 

management,  notably:  what  are  the  biological  conditions  necessary  for  sustainability, 

what  are  the  costs  of  achieving  this  goal,  and  what  are  appropriate  regulatory  mechanisms 

for  inducing  sustainable  behavior.  For  reasons  of  tractability,  numerous  assumptions  are 

made  in  developing  the  biological  and  economic  models.  Consequently,  although  the 

models  are  derived  in  rigorous  empirical  analysis,  the  simulation  results  are  indicative, 

not  definitive.  It  is  hoped,  nonetheless,  that  the  results  of  this  study  will  inform 

management  and  policy  formulation  for  the  Chiquitano  forests  and  enrich  important 

research  concerning  the  biological  and  economic  requirements  for  the  sustainable 

management  of  tropical  forests. 


CHAPTER 2 

MULTINOMIAL  LOGIT  ESTIMATION  OF  A  MATRIX  GROWTH  MODEL  FOR 

TROPICAL  DRY  FORESTS  OF  EASTERN  BOLIVIA 

Introduction 

The  development  of  growth  and  yield  models  for  tropical  forests  is  often  a  daunting 
task  because  of  the  complexity  of  these  rich,  diverse  ecosystems,  and  the  paucity  of  data 
enabling  robust  model  estimation.  The  difficulties  of  representing  the  dynamics  of 
species  interaction,  recruitment,  and  response  to  disturbance  are  compounded  when  the 
model  is  intended  to  serve  goals  of  management  planning,  in  addition  to  scientific 
simulation.  In  the  present  study,  I  develop  a  matrix  model  of  forest  growth  and  yield  for 
tropical  dry  forests  of  eastern  Bolivia,  the  product  of  several  years  of  tree  growth 
measurements  conducted  by  members  of  the  Bolivian  forest  management  project 
(BOLFOR).  The  model  was  developed  to  respond  to  the  need  for  meticulous  estimation 
of  timber  returns  and  forest  impacts  from  polycyclic  harvesting.  The  study  contributes  to 
forest  growth  modeling  science  through  its  novel  use  of  multinomial  logistic  regression 
methods  to  estimate  a  five-guild  matrix  growth  model  for  a  Chiquitano  tropical  dry 
forest. 

Matrix  models  were  first  formulated  by  Leslie  (1945)  for  animal  populations  and 
later  modified  by  Usher  (1966,  1969)  for  managed  forests.  Buongiorno  and  Michie 
(1980)  extended  this  work  by  introducing  density-dependent  recruitment  to  solve  the 


problem  of  exponential  growth  of  the  tree  population  in  the  Leslie  and  Usher  models. 
Their  model  formulation  provided  the  basis  for  a  large  number  of  later  studies  by 
Buongiorno  and  colleagues  that  integrated  species  stratification,  forest  regulation,  and 
economic  and  ecological  criteria  with  economic  optimization  models  founded  upon 
density-dependent  matrix  growth  model  estimation  (e.g.,  Buongiomo  and  Lu  1990,  Lu 
and  Buongiorno  1993,  Buongiorno  et  al.  1994,  Ingram  and  Buongiorno  1996,  Boscolo 
and  Buongiorno  1997). 

Matrix  models  are  based  on  projections  of  whole  stands,  as  opposed  to  individual 
trees.  Trees  are  commonly  aggregated  according  to  their  size  class  and  perhaps  species 
group  so  that  the  entire  stand  may  be  represented  by  a  vector  (or  several  vectors)  of  tree 
population  per  class  at  a  given  time.  Matrix  models  are  based  on  a  system  of  linear 
difference  equations,  describing  the  change  in  tree  populations  per  size  class  over  discrete 
time  periods.  Accordingly,  for  any  given  population  value  in  a  size  class  (y,)  at  time 
I  =  0,  where  the  function/describes  the  change  of  y,  over  a  discrete  time  period  for  all 
periods  {y„  =fit,  y«.;),  1  =  1 ,  2,  3 . . .  T) ,  there  exists  a  uniquely  determined  function  yu  that 
is  a  solution  of  the  equation  and  has  the  given  value  for  /  =  0  (Sydsaeter  and  Hammond 
1995). 

An  important  advantage  of  the  matrix  formulation  is  its  relative  simplicity  and  ease 
of  integration  into  linear  programming  optimization  models  (Buongiorno  and  Michie 
1980).  Classifying  the  growing  stock  relative  to  only  a  few  parameters  greatly  simplifies 
the  modeling  of  complex  uneven-aged,  multi-species  stands.  The  linear  form  of  the 
model  and  stand  classification  into  size/species  groups  with  consistent  behavior  greatly 


simplifies  model  calibration  and  inferences  concerning  forest  evolution  (Lu  and 
Buongiorno  1993). 

Transition  probabilities  are  commonly  estimated  by  simple  mean  proportions  of 
movement  calculated  from  observations  of  the  animal  or  tree  populations  (Leslie  1945, 
Usher  1966,  Buongiorno  and  Michie  1980).  Logistic  regression  to  predict  the  probability 
that  an  individual  tree  will  be  'successful'  in  an  event  with  a  binary  outcome  has  been 
used  to  predict  mortality  (Hamilton  and  Edwards  1976,  Buchman  et  al.  1983), 
regeneration  (Johnson  1984,  Ferguson  et  al.  1986),  and  diameter  increment  (Lowell  and 
Mitchell  1987,  Vanclay  1991).  The  logistic  function  constrains  probability  predictions  to 
the  interval  (0,1),  provides  a  binomial  distribution  of  errors,  and  allows  the  use  of 
individual  tree  data  rather  than  class  means  for  tree  size  and  other  explanatory  variables 
(Vanclay  1994).  Moreover,  the  logistic  equation  is  robust  in  the  presence  of  outliers  and 
decrements  in  the  data  (Vanclay  1991). 

In  the  present  study,  I  employ  multinomial  logistic  (MNL)  regression  to  estimate 
the  transition  probabilities  of  a  matrix  growth  model  for  a  dry  forest  of  the  eastern 
Bolivian  lowlands.  The  MNL  estimates  of  transition  probabilities  are  derived  from  tree 
and  stand  attributes  influencing  growth  and  mortality.  Lowell  and  Mitchell  (1987) 
developed  a  logistic  model  that  estimated  growth  and  mortality  simultaneously,  though 
these  were  portrayed  as  binary  choices;  specific  growth  proportions  were  applied  to  the 
"successful"  population  in  a  subsequent  model  step.  Rather  than  estimating  the 
probability  of  "success"  for  a  binary  outcome,  such  as  mortality  or  survival,  the 
multinomial  formulation  allows  estimation  of  the  probability  of  one  of  three  events 
occurring  for  a  size-class  and  species  group  of  trees  during  a  given  growth  interval: 


•  Mortality 

•  Remaining  in  the  same  size  class  or  "stability" 

•  Moving  up  one  size  class  or  "upgrowth" 

This  approach  differs  from  previous  work  in  its  use  of  tree  and  forest  characteristics 
to  simultaneously  estimate  the  probability  of  these  three  outcomes.  The  approach  is 
advantageous  in  allowing  deterministic,  stochastic,  and  dynamic  prediction  of  forest 
evolution,  while  preserving  the  simple  linear  form  of  matrix  models  that  permits  their 
straightforward  integration  into  optimization  studies.  Moreover,  it  is  expected  that  the 
transition  probabilities  resulting  from  MNL  estimation  are  more  smoothly  distributed 
across  size  classes  relative  to  those  resulting  from  simple  proportional  estimation  from 
observed  forest  populations,  which  are  likely  to  be  much  more  erratically  distributed 
given  limited  population  samples  and  time  spans  of  observations. 

The  MNL  transition  model  allows  stochastic  simulation  of  transition  probabilities 
while  preserving  the  stationarity  of  a  deterministic  model  form.  The  model  may  thus 
incorporate  the  biological  uncertainty  inherent  in  forest  growth  estimation  and  may  allow 
for  a  more  robust  statistical  comparison  of  results.  Compatible  deterministic  and 
stochastic  forms  permit  more  robust  statistical  analyses,  while  preserving  the  efficiency 
of  a  deterministic  model  (Vanclay  1991). 

An  important  limitation  of  the  transition  matrix  models  developed  by  Leslie  (1945), 
Usher  (1966, 1969),  and  Buongiorno  and  colleagues  (e.g.,  Buongiorno  and  Michie  1980) 
is  the  assumption  of  stationarity—  that  the  transition  probabilities  for  a  given  size/species 
class  remain  constant  over  time.  Estimations  of  steady-state  stocking  and  harvest  levels 
are  feasible;  and  with  stationarity,  a  unique  solution  for  sustainable  harvesting  and 
growing  stock  levels  can  be  found  (Buongiorno  and  Michie  1 980).  The  stationarity 
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assumption  that  upgrowth  and  mortality  probabilities  remain  constant  over  time  may  be 
violated  in  managed  stands,  however  (Johnson  et  al.  1 99 1 ,  Favrichon  1998).  Forest 
succession  is  likely  a  non-stationary  process.  Thus,  probabilities  of  growth  may  be  better 
represented  in  a  functional  rather  than  constant  form  (Usher  1 979).  The  impacts  of 
harvesting  and  of  inter-harvest  silvicultural  treatments  on  established  tree  growth 
therefore  cannot  be  simulated,  given  the  stationarity  assumption.  Stationarity  may  not 
well  describe  actual  stand  behavior  if,  for  instance,  the  growth  of  larger  trees  responds  to 
structural  changes  caused  by  harvesting,  as  observed  in  humid  tropical  forests  of 
Suriname  and  Brazil  (Graaf  1986,  Jonkers  1987,  Silva  1989). 

The  MNL  model  estimation  enables  dynamic  projections  of  forest  growth  by 
reestimating  transition  probabilities  at  each  time  interval,  thus  overcoming  the 
stationarity  limitation  of  traditional  matrix  models.  The  MNL  model  thus  allows  for 
simulation  of  forest  development  dynamics,  with  transitions  influenced  by  changing 
stand  characteristics.  Decisions  concerning  harvesting  and  other  silvicultural  treatments, 
as  well  as  the  evolution  of  stand  structure  and  density,  may  be  considered  in  iterative 
projections  of  stand  development  at  different  time  intervals.  Stochastic  estimation  may 
also  be  employed  in  the  dynamic  model  form. 

Study  Site 

The  matrix  model  is  estimated  from  permanent  sample  plot  (PSP)  data  collected  in 
two  400  hectare  (ha)  forest  blocks  near  the  community  of  Las  Trancas  (16°1 3'S, 
61°50'W)  in  the  Lomerio  region  south  of  Conception,  Santa  Cruz,  Bolivia.  The 
seasonally  dry  tropical  forests  of  Lomerio  lie  in  a  transition  zone  between  the  humid 
forests  of  the  southern  rim  of  the  Amazon  basin  and  the  thorn  scrub  of  the  Gran  Chaco  on 
the  southwestern  edge  of  the  Brazilian  shield  (Killeen  et  al.  1998).  The  forests  are  typical 


of  the  Chiquitano  dry  forests  of  eastern  Bolivia.  Seasonal  mean  temperature  in  the  region 
is  24.3°C  with  a  mean  annual  precipitation  of  about  1 100  mm  and  an  acute  dry  season 
from  May  to  October.  Soils  of  the  region  are  primarily  Inceptisols  (shallow)  and  Oxisols 
(deep),  found  in  four  distinct  mapping  units:  hilltops,  upper  slopes,  lower  slopes,  and 
valley  bottoms  (Iporre  1996).  The  natural  vegetation  is  classified  as  tropical  dry  forest 
(Holdridge  1967).  The  undulating  topography  is  dominated  by  low  hills  composed  of 
granite,  gneiss,  and  metamorphic  rocks  of  Precambrian  origin  (Geobold  1981).  Elevation 
varies  between  400  m  and  600  m  asl. 

Three  distinct  habitats  constitute  the  Lomerio  forests:  upland  forest;  riparian  or 
valley-bottom  forest;  and  granitic  outcrops  (inselbergs).  Mature  upland  and  riparian 
forest  canopies  range  from  12  -  18  m  tall  and  are  dominated  by  trees  of  the  Leguminosae, 
Bignoniaceae,  Anacardiaceae,  and  Bombacaceae  (Killeen  et  al.  1 998).  Understory  trees 
are  mostly  Sapindaceae  and  Myrtaceae  (Kennard  et  al.  2002).  The  spiny  ground 
bromeliad  Pseudananas  sagenarius  is  distributed  over  approximately  80%  of  the  forest 
(MacDonald  et  al.  1998)  and  regenerates  vigorously  in  upland  forest  clearings.  Light 
selective  logging  for  Spanish  cedar  (Cedrelafissilis),  Spanish  oak  (Amburana  cearensis), 
and  "morado"  (Machaerium  scleroxylori)  occurred  at  undetermined  intervals  in  the  past 
(Fredericksen  et  al.  2001). 

Methods 

BOLFOR  researchers  arranged  180  PSPs  of  20  x  50  m  in  a  stratified  random  design 
in  the  two  400  ha  forest  blocks  in  1994  (LT94,  Claros  and  Licona  1995)  and  1995  (LT95, 
Killeen  et  al.  1 998).  The  PSP  cover  a  total  forest  area  of  1 8  ha:  8  ha  in  plot  LT94  and  1 0 
ha  in  LT95.  Measurements  of  diameter  at  breast  height  (dbh),  crown  and  stem  quality, 
crown  position,  liana  infestation,  damage,  and  mortality  were  conducted  bi-annually. 
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Final  PSP  measurements  for  this  study  were  conducted  in  July  2001 .  The  forest  blocks 
were  lightly  logged  1  year  after  their  installation.  Mean  harvesting  intensities  in  the  PSPs 
were  2.72  trees  per  ha  or  a  basal  area  of  approximately  0.68  m2/ha. 

The  density  of  trees  >  10  cm  dbh  was  estimated  at  418  stems/ha  in  the  initial  PSP 
inventories  of  LT94  and  437  stems/ha  in  LT95,  indicating  that  the  Lomerio  forest  canopy 
is  relatively  open  (Killeen  et  al.  1998).  Of  6005  trees  >  10  cm  dbh  initially  inventoried 
and  monitored  for  growth  in  the  PSPs,  5273  (87.8%)  were  retained  for  growth  modeling. 
The  remaining  trees  were  purged  from  the  data  set  because  of  irregular  stem  form  (8.2%), 
apparent  measurement  error  (3.6%),  or  unknown  species  (0.4%).  Trees  of  irregular  form 
included  species  of  Arecaceae  and  Cactaceae,  as  well  as  species  that  change  stem  form 
with  growth  (Acosmiun  cardenasii,  Aspidosperma  spp.,  Chorisia  speciosa,  and  Ficus 
gomelleira).  The  discarded  trees  were  included  in  estimates  of  plot  BA  employed  in 
model  regressions.  Growth  increments  were  recorded  for  trees  10  to  122  cm  dbh. 
Adjusted  for  season  of  measurement,  the  PSP  growth  data  cover  82  months  (LT94)  and 
68  months  (LT95).  Annual  growth  rates  ranged  from  0.00-1 .83  cm/y. 

Given  the  heterogeneity  of  species  attributes  and  of  forest  habitats,  a  grouping  of 
tree  species  into  ecological  guilds  was  conducted  to  enable  more  precise  modeling  of 
growth,  recruitment,  and  mortality  (Appendix  A).  Ninety-three  tree  species  were 
grouped  into  five  guilds  defined  by  forest  habitat  and  shade  tolerance  of  regeneration 
(Pinard  et  al.  1999).  Shade  tolerance  of  regeneration  is  expected  to  be  a  robust  measure 
of  differences  between  guilds  (Grubb  1977).  For  cases  in  which  the  shade  tolerance  of 
regeneration  was  not  classified,  mature  tree  shade  tolerance  was  employed  in  guild 
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grouping.  The  five  guilds  retained  for  development  of  growth  and  yield  models  for 

Lomerio  are  defined  as  follows: 

Guild  1:  riparian  species,  shade  tolerant 

Guild  2:  riparian  species,  shade  intolerant 

Guild  3:  early  successional  and  granitic  outcrop  species,  shade  intolerant 

Guild  4:  mature  upland  forest  canopy  species  and  generalists,  shade  intolerant 

Guild  5:  mature  upland  canopy  species  and  generalists,  shade  tolerant 


The  growth  model  for  Lomerio  forests  is  founded  upon  a  characterization  of  the 
tree  population  by  structure  and  composition.  Classification  of  trees  by  diameter  class  is 
utilized  to  describe  forest  structure  and  a  grouping  of  species  by  ecological  guild  is 
conducted  both  to  portray  floristic  composition  and  to  enable  more  precise  modeling  of 
regeneration  and  disturbance  events.  Diameter  classes  are  defined  in  five-centimeter 
intervals  from  10  cm  to  85+  cm  dbh  and  the  growth  interval  9  is  defined  as  a  five-year 
period.  The  growth  model  is  specified  as 


y<+a  =  G(y,  -  h,)  +  r, 
where 


(2-1) 


<:, 


yi=\yu,y2i,  — >y«J' 

h,=  [h„,  h2/,  ...,  hm,]'  and 
ft™  [r\h  Tn,  ...,  rm,]' 
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Trees  in  the  stand  are  divided  into  a  finite  number  of  diameter  classes  (n)  and 
guilds  (m)  with  yljt  the  number  of  trees  prior  to  harvest  in  the  i*  diameter  class  of  guild; 
at  time  t.  The  total  density  of  guild;  at  time  t  is  defined  as  the  column  vector 

yjt  =  \s»  ]■ 

The  number  of  trees  harvested  per  diameter  class  (I)  in  a  given  guild  (/)  is  described 
by  the  column  vector  h,,  =  [h,,,]. 

The  probabilities  of  upgrowth  or  movement  from  one  diameter  class  to  the  next 
during  a  5-year  growth  interval  (0)  are  expressed  as  a  matrix  (G)  of  transition  sub- 
matrices,  and  used  to  predict  change  for  a  specific  time  interval  (y0). 


The  matrix  G,  is  composed  of  transition  probabilities  for  trees  of  guild;',  which 

define  the  movement  of  trees  into  size  classes  (i)  during  the  period  0  as  follows: 

a,j         is  the  probability  that  a  tree  in  size  class  I  will  remain  alive  and  in  size 
class  /,  or  stability; 

b,j        is  the  probability  that  a  tree  in  size  class  t-1  will  remain  alive  and  grow 
into  i  from  size  class  ;'-l ,  or  upgrowth; 

Cj,         is  the  probability  that  a  tree  in  size  class  i  will  die  during  the  period,  or 
mortality. 

The  probabilities  are  related  by  the  following  equations: 
%  +  b,y  +  c  ij  =1     for  /  <  n 


,+Cs-l 


for  i m  n 
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The  MNL  growth  model  estimates  the  probability  of  an  individual  tree  observing 

one  of  these  three  transition  alternatives  during  a  five-year  growth  interval.  Mortality 

includes  both  natural  death  and  harvest-induced  death,  which  is  observed  for  those  trees 

immediately  killed  during  harvest  operations  and  those  damaged  during  harvesting  that 

later  die.  Transition  probabilities  are  estimated  as  a  function  of  tree  attributes  and  stand 

characteristics,  described  as  follows: 

DBH  tree  size  class  (dbh  in  5  cm  increments) 

RD/BA      tree  relative  diameter  (tree  dbh/mean  dbh  of  stand)  divided  by  total 
post-harvest  basal  area  per  hectare  (m2/ha) 

TPH  total  number  of  trees  per  hectare  prior  to  harvest 

HBA  harvest  basal  area  per  ha  (m2/ha) 

G  ;  dummy  variable  representing  gui  Id  j,  j  =  2 . . .  5 

Gy*DBH  guild  dummy  variable  multiplied  by  size  class 

Size  class  and  relative  diameter  variables  portray  the  relative  dominance  of  each 
tree  in  its  immediate  forest  environment.  Tree  population,  harvest  basal  area  and  post- 
harvest  basal  area  describe  the  density  and  competition  specific  to  each  plot.  The  dummy 
variables  G,  and  G/DBH  describe  behavioral  differences  among  ecological  guilds, 
attributable  to  their  sensitivity  to  light  and  habitat.  G,  are  intercept  shifters,  while 
G/DBH  change  the  slope  of  the  growth  curve  relative  to  guild  and  size  class  variables. 
The  dummy  variables  have  the  effect  of  increasing  upgrowth  probabilities  for  shade 
intolerant  Guilds  2,  3,  and  4  and  of  decreasing  the  upgrowth  probabilities  for  shade 
intolerant  Guild  5.  The  dummy  variables  adjust  transition  probabilities  for  Guilds  2-5 
relative  to  Guild  1. 
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Transition  probabilities  for  each  guild  are  estimated  according  to  Greene  (2000,  p. 
858).  The  probability  that  tree ;  will  transition  to  state y  is 


Prob(7ree,  =  j)  =  -£ for/  =  0,  1 ,  2;  (2-2) 


where, 

0  =  mortality,  1=  stability,  2  =  upgrowth 

It  is  assumed  for  analytical  convenience  that  the  data  conform  to  the  Independence 
of  Irrelevant  Alternatives  (IIA)  property.  This  property  stipulates  that  "the  ratio  of  the 
probabilities  of  choosing  any  two  alternatives  is  independent  of  the  attributes  of  any  other 
alternative  in  the  choice  set"  (Hausman  and  McFadden  1984,  p.1221).  More  simply 
stated,  IIA  requires  that  the  equation  errors  are  independent  across  states  or  choice 
alternatives. 

The  marginal  effect  of  attribute  x,  on  the  transition  probability  to  state/  is  a 
function  of  the  beta  vector  for  other  states.  Following  Greene  (2000,  p.  861),  the 
marginal  effect  of  x,  is  defined  as 
dP, 


dx,    '   ' 


0j-LW> 


(2-3) 


Recruitment  (the  number  of  live  trees  that  grow  into  the  smallest  diameter  class 
during  a  5-year  interval  t  to  t+S)  is  expressed  as  a  function  of  the  total  tree  population 
after  harvest.  A  negative  relationship  is  expected  between  recruitment  and  residual  tree 
population.  The  recruitment  function  is  of  the  form 

r,=rf+«i  £0V-A„)  (2-4) 
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where  r,  is  the  number  of  recruits  into  the  1 0  cm  size  class  and  d  and  e  are  regression 
parameters  to  be  estimated.  Recruitment  is  estimated  for  the  aggregate  tree  population, 
as  data  did  not  permit  the  estimation  of  statistically  significant  recruitment  models  for 
each  guild.  Recruitment  levels  for  each  guild  are  allocated  proportionately,  according  to 
the  guild  proportion  of  the  total  stand  population  observed  prior  to  harvest.  This  method 
is  imprecise,  but  an  unfortunate  imposition  of  limited  recruitment  data. 

The  population  per  guild  at  time  t+B  is  determined  by  the  situation  at  time  t,  the 
harvest  during  the  growth  interval  (0),  and  the  recruitment  during  this  interval  by 
equations  for  each  of  the  n  guilds 

yU(+8  =  rj,+  a,/yw-h,;,) 

V2/I+9  =  t>2,{y  y,  -  h  y,  )  +  &2j{y2jl  -  kyt) 

y  3/t+s  =  b3/<y2,7-  h2/,)+  a3,{y3,,-  hv) 

Vn/i+e  "  b^(y„-i;/-  W-y,)  +  aw<y„y,-  h^,) 

Growth  may  be  assumed  stationary,  and  long-term  projections  of  stand  growth 
from  initial  conditions  (yo)  made  deterministically  as 

^=G^„+£gV0  (2-5) 

for  y  growth  periods  of  length  6  (Buongiomo  and  Michie  1980).  Alternately,  multiple 
growth  intervals  (yff)  may  be  estimated  iteratively  to  more  closely  reflect  the  dynamics  of 
forest  growth,  with  reestimation  of  the  transition  probabilities  at  each  5-year  interval. 

Results  and  Discussion 
Growth  Model  Estimation 

Multinomial  logit  (MNL)  estimation  resulted  in  a  significant  matrix  model 
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(P  <  0.001)  from  5273  observations  (Table  2-1).  The  marginal  effects  of  tree  and  stand 

attributes  (*,)  on  the  transition  probabilities  are  estimated  at  the  mean  of  Xt  (Table  2-2). 

Table  2-1 .  Maximum  likelihood  estimates  (MLE)  of  transition  model  parameters  *• b 

MLE  Parameter 

Variable 


Stability 


Upgrowth 


DBH 

0.11899 

(0.01126) 

DBH2 

-0.00074 

(0.00016) 

RD/BA 

-9.92015 

(2.58866) 

TPH 

-0.00009* 

(0.00021) 

HBA 

-0.09417 

(0.01693) 

G2 

0.30933* 

(0.77065) 

G3 

0.87904 

(0.40223) 

G4 

0.87814 

(0.27010) 

Os 

0.76966 

(0.19444) 

G2*  DBH 

-0.00239* 

(0.02934) 

G3*  DBH 

-0.03752 

(0.01509) 

G4*  DBH 

-0.05345 

(0.01047) 

G5*DBH 

-0.03654 

(0.00957) 

0.05704 
(0.01285) 

-0.00011* 
(0.00017) 

-8.03101 

(2.99525) 

-0.00096 
(0.00025) 

-0.06784 
(0.01966) 

0.75517* 
(0.82540) 

0.77800f 
(0.46470) 

1.89981 
(0.28991) 

0.88485 
(0.22188) 

-0.00334* 
(0.03114) 

-0.03436 
(0.01660) 

-0.06214 
(0.01137) 

-0.03347 
(0.01072) 


X^^d-f.) 
lnl 


281.3135 
-4590.429 


Notes:     a.  Asymptotic  standard  errors  in  parentheses 

b.  Parameters  for  mortality  are  set  to  zero  in  MNL  estimation 
*  Significant  to  P<  0.10 
Not  significant 
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Table  2-2.  Marginal  effects  of  attributes  x,  on  transition  probabilities  j ' 


Marginal  Effect 


Variable 

Mortality 

Stability 

Upgrowth 

DBH 

-0.01205 
(0.00123) 

0.01911 
(0.00181) 

-0.00706 
(0.00155) 

DBH2 

0.00007 

(0.00002) 

-0.00015 
(0.00002) 

0.00009 

(0.00002) 

RD/BA 

1.10026 
(0.28906) 

-1.13126 
(0.45397) 

0.03100* 

(0.39240) 

TPH 

0.00004* 
(0.00002) 

0.00012 
(0.00004) 

-0.00015 
(0.00003) 

HBA 

0.01020 
(0.00188) 

-0.01192 
(0.00315) 

0.00173* 
(0.00273) 

G2 

-0.04911* 
(0.08785) 

-0.03591* 
(0.10680) 

0.08501* 
(0.08313) 

G3 

-0.09944 
(0.04540) 

0.09088* 
(0.06616) 

0.00856* 
(0.05712) 

G4 

-0.13225 
(0.03000) 

-0.06753* 
(0.04113) 

0.19978 
(0.03297) 

Qi 

-0.09304 
(0.02184) 

0.05085* 
(0.03068) 

0.04218* 
(0.02626) 

G2*  DBH 

0.00031* 
(0.00338) 

-0.0007* 
(0.00369) 

-0.00023* 
(0.00275) 

G3*  DBH 

0.00428 
(0.00171) 

-0.00372* 
(0.00220) 

-0.00056* 
(0.00179) 

G4*  DBH 

0.00648 
(0.00117) 

-0.00343 
(0.00145) 

-0.00305 
(0.00116) 

G5*DBH 

0.00417 
(0.00108) 

-0.00362 
(0.00133) 

-0.00055* 
(0.00111) 

Notes:     a.  Asymptotic  standard  errors  in  parentheses 
+  Significant  to  P<  0.10 
Not  significant 

The  MNL  estimated  transition  probabilities  are  not  significantly  different 

(P  <  0.05)  from  the  observed  probabilities  derived  from  Lomerio  PSP  data  by  calculation 

of  simple  mean  proportions  of  movement  per  size  class  and  guild  (Table  2-3). 
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As  an  example,  the  resulting  MNL  transition  matrix  (G5MNL)  estimated  for  stand 
conditions  after  the  observed  harvest  and  the  transition  matrix  generated  by  simple 
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proportional  estimation  (G50)  for  Guild  5  with  a  five-year  interval  and  5  cm  dbh  classes  is 
described  below.  Guild  5  is  selected  for  comparison  of  observed  and  predicted  transition 
matrices  because  it  had  the  highest  number  of  observations  (2762). 

As  hypothesized,  the  MNL  transition  matrix  results  in  a  smoother  distribution  of 
transition  probabilities  across  size  classes  relative  to  that  generated  by  simple 
proportional  estimation  (Figure  2-1).  The  MNL  estimation  corrects  for  variance  in  the 
forest  data  and  for  sample  error  effects  on  matrix  elements  (Usher  1976),  such  as  the 
elevated  probabilities  of  movement  in  large  size  classes  generated  by  proportional 
estimation  from  a  very  small  sample  of  large  trees. 

1.20 
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Figure  2-1 .  Observed  and  predicted  stability  probabilities  per  dbh  class  for  Guild  5 

Recruitment 

The  OLS  estimated  recruitment  function  (r)  observes  the  expected  negative 
relationship  between  recruitment  and  total  trees  per  hectare  (TPH)  (Table  2-4).  Although 
the  fit  is  quite  poor,  the  variables  and  the  model  are  all  significant  to  P  <  0.05. 
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Table  2-4.  Ordinary  least  squares  (OLS)  estimate  of  recruitment  model 

Variable  Parameter  Std.  Error 


Intercept 

152.775 

27.213 

TPH 

-0.117 

0.054 

AdjR2 

0.021 

— 

Ftest 

4.583 

— 

Static,  Deterministic  Modeling 

In  static,  deterministic  form,  transition  probabilities  estimated  per  guild  at  the  initial 
stand  state  are  maintained  throughout  the  projection  horizon.  The  MNL  model  of  forest 
evolution  after  harvest  shows  an  increase  in  total  tree  population  and  a  leveling  of  the 
population  distribution  across  size  classes  with  increasing  age  (Table  2-5). 
Table  2-5.  Deterministic  prediction  of  forest  evolution,  aggregate  population 


DBH  Class 

Post-harvest 
(trees/ha) 

Projected  (trees/ha) 

(cm) 

year  25 

year  50 

year  100 

10 

170.2 

191.3 

189.6 

187.7 

15 

141.2 

129.4 

129.3 

127.7 

20 

82.7 

88.6 

90.0 

88.8 

25 

47.7 

57.9 

61.9 

62.0 

30 

33.3 

36.1 

40.5 

42.8 

35 

26.5 

23.2 

25.3 

28.7 

40 

19.7 

16.2 

15.8 

18.3 

45 

11.8 

11.2 

10.2 

11.2 

50 

8.6 

7.3 

6.6 

6.8 

55 

4.6 

4.3 

4.1 

4.1 

60 

3.3 

2.6 

2.5 

2.5 

65 

1.9 

1.6 

1.5 

1.6 

70 

1.3 

1.0 

0.9 

1.0 

75 

0.5 

0.6 

0.6 

0.6 

80 

0.7 

0.3 

0.4 

0.4 

85 

0.9 

0.3 

0.3 

0.3 

Total 

554.8 

571.7 

579.3 

584.6 

Comparison  of  model  projections  using  the  observed  transition  probabilities  (G  ) 
and  the  MNL  estimated  probabilities  (GMNL)  reveals  important  differences  in  the 
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resulting  stand  distribution  estimates  (Figure  2-2).  Whereas  the  MNL  estimated  matrix 
model  predicts  a  leveling  population  distribution  with  increasing  age,  the  matrix  model 
derived  from  observed  transition  probabilities  predicts  a  distribution  similar  to  that 
observed  after  the  initial  harvest  and  a  greater  decline  in  stems  of  merchantable  size  (>  40 
cm  dbh).  The  MNL  model  predicts  a  greater  number  of  trees  in  the  merchantable  20  to 
40  cm  dbh  classes  relative  to  the  observed  prediction  due  to  the  "smoothing"  effect  of  the 
regression  on  transition  probabilities  across  size  classes.  While  the  observed  matrix 
model  exhibits  considerable  variance  in  transition  probabilities,  especially  in  those  guilds 
with  limited  samples  (Guilds  2  and  3),  the  MNL  model  corrects  for  sample  variance.  The 
MNL  model  predicts  greater  stability  of  stems  in  the  20  to  40  cm  range  and  greater 
upgrowth  of  smaller  stems  into  these  size  classes  relative  to  the  observed  model,  with  less 
variance  of  these  transition  probabilities  among  similar  size  classes. 
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Figure  2-2.  Static,  deterministic  prediction  of  forest  evolution  with  observed  (G  )  and 
MNL  estimated  (G     L)  transition  matrices 
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Static,  Stochastic  Modeling 

The  MNL  estimation  of  transition  probabilities  permits  stochastic  simulation  of 
stand  evolution,  drawing  upon  the  standard  errors  of  regression  parameters.  The 
stochastic  approach  considers  biological  uncertainty  in  projections  of  stand  growth, 
providing  a  confidence  interval  of  expected  population  values  for  a  future  stand  state. 
Stochastic  simulation  has  the  benefit  of  enabling  more  rigorous  statistical  comparison  of 
results  from  alternative  simulations  and  permitting  sensitivity  analyses.  The  stochastic 
element  can  be  introduced  into  static  projections,  for  which  the  initial  transition 
probability  estimates  are  maintained  throughout  the  projection  horizon  (Figure  2-3)  or  in 
dynamic,  iterative  projections. 
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10     15     20    25     30     35     40    45     50     55     60     65     70    75     80     85 


DBH  class  (cm) 


Figure  2-3.  Stochastic,  stationary  estimation  of  the  forest  population  (all  guilds)  25  years 
after  the  observed  harvest  (SD  =  standard  deviation,  perc  =  percentile) 

Dynamic  Modeling 

Changes  in  the  growth  dynamics  of  the  forest  over  time  may  be  considered  by 
iterative  estimation  of  forest  evolution  with  the  re-estimation  of  transition  probabilities  at 
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each  5-year  growth  interval.  Forest  response  and  recovery  trajectories  following 
disturbances  caused  in  harvesting  or  silviculture  may  be  more  precisely  estimated  in 
dynamic  form.  As  an  example,  estimated  stand  populations  and  their  distribution  twenty- 
five  years  after  a  simulated  harvest  of  40%  of  all  merchantable  stems  (>  40  cm  dbh; 
Guilds  1, 2, 4,and  5)  are  quite  distinct  in  static  and  dynamic  model  forms  (Figure  2-4). 
The  25-year  growth  horizon  simulated  in  this  example  was  selected  merely  for 
demonstration  purposes.  Changes  in  the  stand  population,  basal  area,  and  harvest 
variables  over  time  result  in  modest  changes  in  the  transition  probability  estimates 
generated  by  the  MNL  model,  which  may  more  precisely  reflect  forest  dynamics  and 
disturbance  effects. 
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Figure  2-4.  Comparison  of  static  and  dynamic  model  estimates  of  tree  population 
distribution  25  years  after  harvesting  40%  of  merchantable  stems 

That  the  dynamic  model  predicts  more  stems  in  10-40  cm  dbh  classes  relative  to 

static  model  (Figure  2-4)  reveals  an  important  advantage  of  the  dynamic  prediction  form. 

In  the  static  model,  the  transition  probabilities  estimated  for  the  initial  harvest  are 
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maintained  throughout  the  25-year  projection  period.  Consequently,  the  high  mortality 
probabilities  estimated  after  a  40%  proportional  harvest  are  retained  at  each  5-year 
growth  interval  in  the  static  model,  reducing  the  residual  population  at  each  interval.  The 
dynamic  model  avoids  this  modeling  limitation  by  recalculating  transition  probabilities  at 
each  growth  interval,  according  to  the  state  of  the  stand  at  these  5-year  intervals. 
Accordingly,  the  high  mortality  probabilities  estimated  immediately  after  the  simulated 
40%  harvest,  are  not  maintained  throughout  the  projection.  Instead,  mortality 
probabilities  decline  with  increasing  time  after  this  harvest,  as  the  stand  returns  to 
relatively  stable  conditions  compared  to  the  disturbed  state  created  by  harvesting. 
An  Application 

The  benefits  of  the  MNL  model  form  relative  to  a  matrix  model  generated  from 
observed  transition  probabilities  are  best  revealed  in  model  application.  Following 
Buongiorno  and  Michie  (1980),  the  undisturbed,  steady-state  or  "climax"  equilibrium 
(ya)  may  be  calculated  for  a  matrix  growth  model  as 

y°=(I-G)-V  (2-6) 

Expectations  for  the  undisturbed  steady-state  are  that  the  tree  population  converge  to  a 
negative  exponential  distribution,  typical  of  uneven-aged,  mixed  species  forests  and  that 
with  greater  length  between  disturbances  this  distribution  gradually  even  out  across  size 
classes.  The  MNL  steady-state  demonstrates  precisely  this  result  (Fig.  2-5).  The  steady- 
state  generated  by  the  matrix  model  derived  from  observed  transition  probabilities 
suggests  a  relative  decline  of  large  stems  with  increasing  age.  Moreover,  an  anomalous 
"bump"  in  the  population  distribution  for  trees  of  size  class  70  cm  dbh  results  from  model 
estimation  errors  imposed  by  the  limited  sample  of  trees  in  large  stem  classes. 


26 


Observed  matrix  model 

160- 

2 

"2  120 " 

6 

I      80  J 

\\ 

40- 
0  - 

1 1 1 ■ 1 1                          1                 " 1 1 

10 


20 


30 


40 


50 


60  70  80 

Dbh  class  (cm) 


Figure  2-5.  Undisturbed,  steady-state  equilibrium  condition  (ya)  estimated  with  observed 
and  MNL-estimated  transition  probability  matrices 

Conclusions 

Multinomial  logistic  (MNL)  estimation  of  the  transition  probabilities  in  matrix 
models  of  forest  growth  offers  important  benefits  relative  to  matrices  derived  from 
proportional,  static  estimates  derived  from  forest  plot  observation  means.  The  MNL 
model  results  in  a  smoother  distribution  of  transition  probabilities  across  size  classes, 
correcting  for  variance  in  the  data  and  model  estimation  errors  imposed  by  limited 
samples.  The  MNL  form  also  allows  for  deterministic  and  stochastic  projection  in  both 
static  and  dynamic  model  forms.  Stochastic  simulation  enables  more  rigorous  testing  of 
hypotheses  concerning  stand  growth  by  accounting  for  fundamental  biological 
uncertainties.  The  static  model  form  allows  greater  ease  of  integration  in  optimization 
studies,  while  dynamic  forms  enable  more  meticulous  simulation  of  stand  evolution 
dynamics  and  the  effects  of  disturbance  events. 


CHAPTER  3 
OPTIMAL  MANAGEMENT  OF  A  CHIQUITANO  TROPICAL  DRY  FOREST 

Introduction 

Sustainable  forest  management  has  drawn  great  attention  in  the  debate  over 

appropriate  strategies  to  achieve  often  competing,  but  necessarily  compatible  objectives 

of  environmental  conservation  and  sustainable  economic  development  in  developing 

nations  with  large  tropical  forest  reserves  (Kishor  and  Constantino  1994,  Dickinson  et  al. 

1996,  Rice  et  al.  1997,  Frumhoff  and  Losos  1998,  Rice  et  al.  2001,  Pearce  et  al.  2003). 

Despite  conclusions  concerning  the  viability  or  appropriateness  of  sustainable 

management  initiatives  within  this  polemic  debate,  the  conditions  necessary  for 

sustainable  management  and,  consequently,  rigorous  analyses  of  the  costs  and  benefits  of 

these  strategies  relative  to  alternative  tropical  forestry  options  have  been  conducted  in 

only  a  few  noteworthy  studies  (Howard  et  al.  1996,  Howard  and  Valerio  1996,  Ingram 

and  Buongiomo  1996,  Boscolo  and  Buongiorno  1997,  Boscolo  and  Vincent  2000).  This 

deficiency  is  primarily  attributable  to  the  paucity  of  statistically  robust  models  of  the 

growth  and  yield  of  mixed-species  tropical  forests.  This  paucity,  in  turn,  is  due  to  the 

common  inadequacy  of  tropical  forest  growth  data,  particularly  that  offering  evidence  of 

forest  responses  to  harvest  interventions.  Furthermore,  a  precise  definition  of  what 

constitutes  SFM  remains  elusive,  but  commonly  sustainability  is  defined  by  one  of  two 

management  objectives: 

•        Maintaining  timber  yield  over  an  indefinite  harvest  horizon  or  sustained-yield 
timber  management  (STY) 
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•        Maintaining  non-declining  timber  yields;  and  sustaining  other  forest  benefits  (non- 
timber  products,  ecological  services,  and  biodiversity) 

The  latter  may  be  more  strictly  defined  as  SFM;  a  possible  distinction  being  that  STY 

may  generate  a  less  sustainable  flow  of  non-timber  products  relative  to  SFM  (Pearce  et 

al.  2003). 

Financial  analyses  comparing  returns  from  more  environmentally  benign,  or 
hypothetically  sustainable,  harvest  prescriptions  with  those  from  conventional  logging 
practices  in  tropical  forests,  reveal  in  some  instances  that  improved  forest  management, 
such  as  reduced-impact  logging,  is  financially  and  ecologically  superior,  though  the 
conditions  required  for  sustainability  are  not  precisely  identified  (Barreto  et  al.  1998, 
Armstrong  2000,  Holmes  et  al.  2002).  Other  studies  estimating  the  reduced  harvest 
levels  and  conservation  easements  necessary  for  SFM  in  natural  tropical  forests  conclude 
that  the  opportunity  costs  of  more  sustainable  harvest  regimes  relative  to  conventional 
logging  or  other  land  uses  are  excessively  high  (Kishor  and  Constantino  1993,  Howard 
and  Valerio  1996,  Hout  1999,  Tay  1999,  Healey  et  al.  2000,  Pinard  et  al.  2000).  Still, 
suitably  indicative  estimates  of  the  financial  costs  and  ecological  benefits  of  STY,  let 
alone  SFM,  remain  elusive  due  to  the  inadequacy  of  dependable  growth  models  and  rare 
attempts  to  integrate  these  models  into  an  appropriate  form  for  rigorous  economic 
analysis  of  forest  management  alternatives. 

The  present  study  aims  to  make  a  modest  contribution  to  resolving  this  inadequacy 
by  producing  a  statistically  robust  economic  comparison  of  management  alternatives  for 
forests  representative  of  the  Chiquitano  tropical  dry  forests  of  eastern  Bolivia.  The 
analysis  draws  upon  a  matrix  model  of  forest  growth  and  yield  and  an  economic 
optimization  model,  integrated  to  enable  the  estimation  of  conditions  necessary  for 
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sustained  timber  yield  and  a  comparison  of  this  optimum  with  optimal  harvest  solutions 
for  alternative  forest  management  regimes.  The  optimization  study  concerns  the 
management  by  large  timber  firms  of  Chiquitano  forest  concessions  and  derives  from 
economic  and  forest  data  collected  in  the  Chiquitano  region  and  the  provincial  capital, 
Santa  Cruz,  Bolivia.  Alternative  management  scenarios  examined  in  the  study  include 
unregulated  or  unconstrained  harvest  (U),  management  under  constraints  imposed  by  the 
Bolivian  forestry  law  (BFL),  sustained  timber  yield  management  (STY),  and 
management  constrained  by  objectives  of  sustaining  timber  yield  and  maintaining  the 
structure  and  composition  of  a  theoretical  "climax  equilibrium"  forest  (SFM).  The 
harvest  horizon  for  all  scenarios  is  set  at  40  years  to  allow  for  comparison  of  the  net 
present  value  of  returns.  The  actual  harvest  horizons  for  all  scenarios  are  longer,  those  of 
the  sustainable  scenarios  STY  and  SFM  being  infinite,  however  the  harvest  horizon  was 
limited  to  40  years  due  to  poor  fit  of  the  recruitment  function  in  the  growth  model,  which 
would  significantly  affect  stocking  and  harvest  levels  for  longer  horizons.  Forty  years 
was  selected  since  it  is  the  legal  concession  period  under  the  Bolivian  forestry  law. 
Harvesting  is  selective,  its  selectivity  defined  by  the  value  and  merchantability  of  trees 
and  constraints  on  their  extraction. 

The  SFM  scenario  does  not  derive  from  mathematical  conditions  describing 
sustained  production  of  non-timber  products  or  ecological  benefit  flows,  but  rather  that  of 
attaining  a  residual  forest  condition  approximating  the  structure  and  composition  of  an 
undisturbed  forest,  as  defined  by  the  growth  model.  The  theoretical  "climax"  condition 
(y°)  or  undisturbed  equilibrium  distribution  is  defined  mathematically 

y"=(l-Gj'r  (3-1) 
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where 

•  I  is  the  identify  matrix 

•  G  is  a  matrix  of  transition  probabilities  defining  the  movement  of  trees  from  one 
size  class  to  the  next  during  a  5-year  growth  interval 

•  r  is  a  recruitment  function  defining  the  number  of  trees  entering  the  smallest  size 
class  during  a  5-year  growth  interval. 

These  and  all  other  variables  used  in  the  present  study  are  defined  in  Appendix  D.  The 
climax  condition  is  defined  solely  by  the  growth  model  and  does  not  assume  that  the 
forest  condition  prior  to  harvesting  is  in  this  undisturbed  equilibrium. 

Conserving  forest  vegetative  structure  and  composition  is  assumed  an  appropriate 
proxy  for  the  maintenance  of  non-timber  and  ecological  benefit  flows  (cf,  Terborgh 
1986,  Hunter  1990).  The  SFM  regime  thus  hypothetically  achieves  the  objectives  of 
SFM  defined  above,  though  empirical  evidence  of  the  maintenance  of  non-timber 
benefits  is  not  presented.  Results  of  the  present  study  are  at  best  indicative,  given  the 
assumptions  required  for  modeling;  however,  they  provide  important  insight  into  the 
conditions  required  for  sustainability  and  the  costs  and  benefits  of  sustainable 
prescriptions  relative  to  alternative  harvest  possibilities. 

Methods 
Forest  Growth  Model 

The  optimization  study  employs  a  five-guild  matrix  model  of  forest  growth,  which 

was  estimated  by  multinomial  logit  (MNL)  regression  of  permanent  sample  plot  (PSP) 

data  collected  in  two  400  hectare  (ha)  forest  blocks  near  the  community  of  Las  Trancas 

(16°13'S,  61°50'W)  in  the  Lomerio  region  south  of  Conception  (cf.,  Chapter  2).  The 

1  The  theoretical  climax  condition  should  arguably  be  presented  in  quotations  or  italics,  given  its 
debatable  nature.  For  reading  ease,  however,  such  formatting  will  not  be  employed  throughout  the  text. 
Climax  is  understood  to  refer  to  the  theoretical  forest  condition  defined  mathematically  as  y°=  (I-G)'r. 
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MNL  matrix  growth  model  projects  the  probabilities  of  upgrowth,  stability,  or  mortality 
of  trees  in  each  of  five  ecological  guilds  during  a  5-year  growth  period  and  predicts 
recruitment  of  trees  into  the  smallest  size  class  (10  cm  dbh)  for  each  guild.  The  linear 
form  of  the  MNL  model  and  the  assumption  of  stationarity  for  the  transition  probabilities 
throughout  the  harvest  horizon  enable  its  straightforward  integration  into  the  present 
linear  optimization  study.  The  Las  Trancas  sites  are  typical  of  the  Chiquitano  dry  forests, 
and  thus  provide  a  representative  stand  for  the  examination  of  management  alternatives  in 
Chiquitano  concessions.  Light  selective  logging  for  Spanish  Cedar  (Cedrelaflssilis), 
Spanish  Oak  {Amburana  cearensis),  and  "morado"  (Machaerium  sckroxylori)  occurred  at 
undetermined  intervals  in  the  past  (Fredericksen  et  al.  2001),  so  they  are  not  in  an 
undisturbed,  climax  state. 
Economic  Data 

The  Bolivian  forestry  law  was  revised  in  1996  (Ley  1 700)  to  more  effectively 
promote  sustainable  management  of  the  nations  forest  resources  (Art.  #1).  The  law  set 
aside  areas  for  conservation  and  protection  purposes,  and  permits  logging  only  in  other 
forest  areas.  Logging  concessions  are  allocated  for  a  40-year  duration  upon  formal 
application  to,  and  approval  by,  Bolivia's  forest  service  (Supehntendencia  Forestal). 
The  law  requires  forest  management  plans  for  all  concessions  and  for  forest 
authorizations  in  private  lands.  Concessions  also  can  be  traded  and  inherited  and  are 
renewable  upon  validation  of  the  concessionaire's  observance  of  sustainable  forest 
management  plans. 

Forest  management  costs  are  drawn  from  surveys  of  four  industrial  forestry  firms 
based  in  Santa  Cruz  and  possessing  logging  concessions  in  Chiquitano  dry  forests. 
Structured  surveys  were  conducted  in  November  2000  and  validated  in  follow-up 
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interview  with  each  of  the  firms  in  July  2001 .  All  costs  and  timber  prices  are  reported  in 
2001  $US  and  are  assumed  constant  throughout  the  management  horizon  (Appendix  B). 
Management  costs  FOB  forest  mill  are  classified  as  one  of  three  types: 

•  Variable  costs  (C)  in  $/m3  for  felling,  skidding  and  log  deck  operations  incurred 
relative  to  the  level  of  harvest  intensity 

•  Fixed  costs  (F)  in  $/ha  incurred  regardless  of  harvest  intensity  at  each  cutting-cycle 
entry  for  planning  and  capital  costs  including  depreciation 

•  Annual  costs  (A)  in  $/ha  paid  throughout  the  harvesting  horizon  regardless  of 
harvest  intensity,  which  include  an  area  fee  of  $l/ha  for  the  concession  (patente) 
and  inscription  fees  to  government  agencies  and  trade  associations 

Variable  costs  are  estimated  at  $12.1  l/m\  fixed  costs  at  $62.51/ha  and  annual  costs 
at  $1.55/ha.  Variable  costs  are  estimated  relative  to  the  efficiency  of  felling  and  skidding 
(m  /hr)  and  the  hourly  costs  of  labor  and  machinery  used  in  these  operations.  Felling 
efficiency  of  3.48  m3/hr  is  the  statistical  mode  for  data  collected  in  Chiquitano  forest 
concessions  of  two-person  felling  teams  in  both  planned  and  unplanned  harvesting 
operations  (Cavero  1998,  Menacho  1999)  and  fit  to  a  gamma  distribution.  Skidding 
efficiency  of  3.68  m  /hr  is  the  statistical  mode  for  data  collected  in  studies  of  operations 
using  rubber-tired  skidders  in  three  forest  concessions  in  Chiquitania  (Alarcon  1997, 
Pafifio  1997,  Crespo  1999)  and  fit  to  a  gamma  distribution. 

A  discount  rate  of  1 8.75  %  is  assumed  for  the  analysis.  This  is  the  mean  real 
interest  rate  on  loans  for  2001  reported  by  Bolivia's  Central  Bank  (Banco  Central  de 
Bolivia  2002).  Market  prices  for  merchantable  stems  per  guild  delivered  to  concession- 
based  timber  mills  (FOB  mill)  are  drawn  from  surveys  of  timber  consultants  and  lumber 
yards  in  Santa  Cruz,  and  from  market  surveys  of  Chiquitania  conducted  in  2001  by 
BOLFOR  researchers  (Appendix  C).  Timber  species  demanded  in  Santa  Cruz 
lumberyards  in  2001  are  assumed  to  comprise  the  merchantable  stock.  Market  value  per 
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guild  (fj)  is  estimated  as  the  weighted  mean  price  for  merchantable  stems  of  Guildy 
(Table  3-1).  This  assumption  greatly  simplifies  market  value  calculations  for  multiple- 
species  harvest.  However,  greater  precision  in  calculating  the  market  value  of  timber  was 
necessarily  sacrificed  for  increased  tractability  of  the  optimization  model. 
Guild  populations  comprise  merchantable  and  non-merchantable  stems.  Harvest  is 
necessarily  constrained  by  the  merchantability  of  the  timber  stock,  as  loggers  will  only 
harvest  those  stems  that  are  of  merchantable  species  and  form.  It  is  assumed  that 
harvesting  is  restricted  to  logs  appropriate  for  saw-timber  milling,  i.e.  stems  of  40  cm  dbh 
and  greater.  The  merchantable  proportion  of  stems  per  guild  (<%)  is  estimated  as  the 
product  of  the  proportion  of  merchantable  species  (sj)  per  guild  and  the  proportion  of 
stems  of  merchantable  form  (ff)  in  that  guild,  based  upon  PSP  data  for  the  sample  forest 
blocks  (Table  3-1).  Trees  of  Guild  3,  those  growing  in  inselberg  habitats,  are  not 
merchantable. 
Table  3-1 .  Market  values  and  merchantable  proportions  of  stems  per  guild 


Market  value  (pj) 

Merchantable  proportion 

Guild 

($US/m3) 

Species  (sf) 

Form(#) 

Stems  (a>,) 

1 

41.28 

0.52 

0.78 

0.41 

2 

62.50 

0.90 

0.96 

0.87 

3 

0.00 

0.00 

0.88 

0.00 

4 

57.03 

0.69 

0.87 

0.60 

5 

104.53 

0.26 

0.84 

0.22 

If  the  forest  structure  at  time  t  is  represented  by  the  vector  y,  =  [y;,,  y^, . . . ,  yra,]', 
forest  volume  is  obtained  by  pre-multiplying  y,  by  the  square  matrix 
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Pi 


Y  = 


¥2 


¥„, 


where  \yj  =  [y^],  j  =  l,...,n;  /  »  l,...,m;  is  the  mean  extracted  volume  of  a  tree  in 

diameter  class  i  and  guild/  The  mean  volume  extracted  in  harvesting  is  estimated  as 

85%  of  the  mean  volume  per  dbh  class  and  guild,  assuming  15%  of  merchantable  volume 

lost  in  logging  waste  (Table  3-2).  The  15%  waste  proportion  is  slightly  higher  than  the 

12.3%  reported  in  production  forests  of  Brazil  for  high  stumps,  split  logs,  and  bucking 

waste  (Holmes  et  al.  2002). 

Table  3-2.  Mean  volume  of  stem  extracted  in  harvest  per  size  class  and  guild  (<|/„) 
Mean  extracted  volume  per  stem  (m3/stem) 


Dbh  class        Guild  1 


Guild  2 


Guild  3 


Guild  4 


Guild  5 


40 

0.65 

0.71 

0.71 

0.74 

0.62 

45 

0.88 

0.91 

0.92 

0.97 

0.83 

50 

1.15 

1.14 

1.16 

1.24 

1.08 

55 

1.47 

1.40 

1.44 

1.54 

1.37 

60 

1.85 

1.68 

1.75 

1.89 

1.70 

65 

2.28 

2.00 

2.09 

2.29 

2.09 

70 

2.78 

2.35 

2.48 

2.73 

2.52 

75 

3.34 

2.74 

2.91 

3.22 

3.01 

80 

3.97 

3.16 

3.37 

3.76 

3.56 

85+ 

4.67 

3.61 

3.88 

4.36 

4.17 

The  net  value  (v;,)  of  an  extracted  stem  of  size ;'  in  guildy  in  2001  $US/m3 
(Table  3-3)  is  thus  calculated  as  the  product  of  the  mean  price  per  guild  (pj)  less  variable 
harvesting  costs  (C)  of  $12.1 1/m3  and  the  mean  extracted  volume  per  stem  of  a  given  size 
class  and  guild  (ys),  assuming  no  economies  of  scale  in  harvesting.  The  net  value  of 
stems  in  Guild  3  is  negative  as  species  have  no  market  value,  thus  harvesting  would 
strictly  impose  costs. 
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Table  3-3.  Net  value  (v„)  of  an  extracted  stem  per  size  class  and  guild 

Net  value  ($/stem) 

Dbh  class       Guild  1  Guild  2  Guild  3  Guild  4  Guild  5 


40 

22.33 

41.96 

-10.11 

39.17 

67.76 

45 

30.11 

53.86 

-13.10 

51.25 

90.34 

50 

39.46 

67.44 

-16.54 

65.33 

117.18 

55 

50.52 

82.74 

-20.47 

81.54 

148.63 

(,0 

63.41 

99.82 

-24.89 

100.01 

185.02 

65 

78.28 

118.75 

-29.84 

120.86 

226.70 

70 

95.27 

139.56 

-35.33 

144.21 

274.00 

75 

114.50 

162.33 

-41.39 

170.18 

327.26 

80 

136.13 

187.11 

-48.04 

198.91 

386.82 

85+ 

160.28 

213.95 

-55.30 

230.50 

453.02 

Optimization  Models 

The  optimization  model  prescribes  maximization  of  the  net  present  value  (NPV)  of 
polycyclic  harvesting  (Z)  subject  to  biological,  market,  policy,  and  sustainability 
constraints.  The  objective  function  does  not  consider  the  costs  of  investment  in  the 
growing  stock  described  by  the  economic  stocking  rule  for  selection  harvest  (Duerr  and 
Bond  1952)  or  soil  expectation  value  (SEV)  used  in  similar  optimization  studies 
(Buongiorno  and  Michie  1980,  Boscolo  and  Buongiomo  1997,  Bach  1999,  Kant  1999). 
These  costs  are  opportunity  costs  of  not  harvesting  merchantable  stems  in  the  initial 
harvest  entry  in  order  to  maintain  the  timber  growing  stock.  The  decision  to  exclude  the 
opportunity  costs  of  maintaining  the  growing  stock  from  the  calculation  of  net  returns  is 
based  upon  the  nature  of  the  forest  property.  As  these  forests  are  allocated  in  concessions 
to  private  industry  under  government  ownership,  the  opportunity  costs  of  conserving  the 
residual  stock  are  public  rather  than  private  costs.  These  opportunity  costs  are  revealed 
in  the  difference  between  unconstrained  and  constrained  harvest  cases  described  below. 
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The  management  problems  are  posed  as  linear  programs  and,  as  such,  stationarity 
is  assumed  for  the  transition  probabilities  of  the  matrix  growth  model  to  permit  a  global 
optimum  solution  to  each  problem.  The  transition  probabilities  estimated  for  the  PSP 
following  their  harvest  in  1995  and  1996  are  selected  as  the  stationary  probabilities  in  the 
guild  growth  matrices,  as  these  estimates  are  expected  to  most  closely  approximate 
growth  in  these  forests.  Four  forest  management  alternatives  are  considered  in  the 
optimization,  notably: 

•  Unconstrained  harvesting  at  40-year  intervals  (U) 

•  Harvesting  under  Bolivian  forestry  law  constraints  at  40-year  intervals  (BFL) 

•  Sustained  timber  yield  management  at  5-year  intervals  (STY) 

•  Sustainable  forest  management  at  5-year  intervals  (SFM) 

A  full  mathematical  description  of  each  alternative  is  presented  in  Appendix  E. 

The  cutting-cycle  between  harvesting  under  regimes  U  and  BFL  is  set  at  40-years, 
the  minimum  interval  required  for  profitable  harvest  (NPV,  >  0),  while  that  of  the  STY 
and  SFM  scenarios  is  5  years,  selected  to  maximize  the  NPV  of  returns  from  these  less 
intensive  harvesting  systems.  Financial  returns  from  the  sustainable  scenarios  (STY  and 
SFM)  do  not  include  any  returns  from  regulation  harvests,  which  would  be  undertaken  to 
convert  the  initial  stand  to  steady-state,  sustained-yield  conditions.  The  harvest  horizon 
for  all  scenarios  was  set  at  40  years  to  allow  for  comparison  of  the  NPV  of  returns. 

The  objective  function  for  U,  BFL  and  STY  cases  prescribes  maximization  of  the 
NPV  of  returns  from  polycyclic  harvesting,  Z^y),  calculated  as  the  present  value  of  a 
terminating  series  at  discount  rate  <5,  with  harvest  horizon  T=  40  years  and  cutting-cycle 
length  /  =  yO,  with  period  yO  being  the  number  of  intervals  (y)  of  5-year  growth  periods 
(ff)  in  a  cutting-cycle.  Cutting-cycles  are  40  years  (r  =  40,  y  =  8)  for  U  and  BFL  and  5 
years  for  STY  and  SFM  (f  =  5,  y  =1).  The  objective  function  is  defined  mathematically 
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All  management  alternatives  are  constrained  by  the  forest  growth  model,  which 
defines  growth  and  recruitment  in  the  residual  forest  between  harvest  intervals.  All  cases 
are  also  constrained  by  absolute  harvest  restrictions,  requiring  that  the  harvest  level  not 
exceed  the  growing  stock  (yt  -  h,  >  0)  and  non-negativity  constraints  (yt,  h,  >  0). 

Forest  growth  constraint  for  cases  U  and  BFL  require  that  the  sum  of  the  harvest 
and  growing  stock  at  time  t  equal  the  residual  stocking  at  time  /-/  plus  any  growth  and 
recruitment  in  this  growing  stock  between  these  cutting-cycles  (Eq.  3-3).  This  forest 
growth  constraint  is  defined  as  follows 


r-i 


h,+y,-G'(yl_l-h,.1)  =  ^Gkr  (3-3) 

Merchantability  restrictions  require  that  the  proportion  of  trees  harvested  per  size 

class  ;'  and  guildy  at  time  t  be  equal  or  inferior  to  the  merchantable  stock  (jtDJ)  present  in 

the  stand.  The  constraint  is  expressed 

\r-»yGV+Vs0  (3-4) 

where  mj  is  the  merchantable  proportion  per  guild. 

For  subsequent  cutting-cycles  of  cases  U  and  BFL,  the  merchantable  proportion  is 

expected  to  decrease,  as  merchantable  volume  is  removed  and  upgrowth  during  the 

cutting-cycle  interval  does  not  replenish  the  merchantable  stock  to  initial  levels.  Thus, 

while  equation  3-4  is  applied  as  a  constraint  to  the  initial  harvest  of  U  and  BFL, 

subsequent  merchantable  stock  constraints  are  comprised  of  the  growing  stock  of  residual 

merchantable  stems  remaining  after  harvest  (f,.:)  plus  any  upgrowth  into  this 
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merchantable  stock  during  the  cutting-cycle  interval.  This  constraint  is  described 
mathematically 

^-G^V.-I^GV^0  (3-5) 

In  addition,  merchantability  standards  require  that  trees  equal  or  exceed  40cm  dbh,  or 

IX  =0  (3-6) 

1-0 

The  Bolivian  forestry  law  constrained  management  alternative  (BFL)  employs  the 
same  objective  function  (Eq.  3-2)  and  is  constrained  by  the  same  forest  growth  (Eq.  3-3) 
and  merchantability  (Eq.  3-4,  3-5)  restrictions  as  the  unconstrained  alternative  (U).  The 
BFL  case  is  constrained  by  the  technical  standards  established  by  the  government  of 
Bolivia  for  all  concessions  and  private  forests  >  200  ha  (Ministerio  de  Desarollo 
Sostenible  y  Planificacion  1998).  The  standards  stipulate  a  minimum  20-year  cutting- 
cycle,  a  maximum  80%  allowable  cut  of  the  merchantable  stock  per  species,  and  a 
minimum  dbh  limit  for  harvested  stems.  The  20-year  cutting-cycle  restriction  is  not 
binding  for  the  present  analysis,  given  that  the  minimum  cutting-cycle  for  positive  NPV 
was  estimated  at  40  years  and  this  interval  assumed  for  the  BFL  scenario. 

The  minimum  dbh  limit  standard  is  not  uniform  for  all  species  of  the  Chiquitano 
dry  forest.  Stems  that  are  more  valuable  generally  have  higher  dbh  limits,  such  as 
Cedrela  fissilis  (60  cm)  and  Amburana  cearensis  (50  cm),  while  the  absolute  minimum 
for  other  stems  is  40  cm  dbh.  The  guild  grouping  employed  in  this  present  model  does 
not  allow  for  precise  application  of  this  diversity  of  dbh  limits.  Rather,  a  minimum  dbh 
of  45  cm  is  assumed  for  all  merchantable  stems  as  a  BFL  constraint  (Eq.  3-7).  Similarly, 
the  80%  allowable  cut  restriction  is  applied  per  guild  and  not  per  species  (Eq.  3-8, 3-9). 
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These  assumptions  decrease  the  precision  of  model  projections,  but  are  necessary  for 
tractability. 

40 

I>,„=0  (3-7) 

1*0 

f>„-0.8f>A,<0  (3-8) 

i=45  /=45 

85  85  85    y-\ 

£%  -0.8XGY,-,  -0.8XZ^G*Vl  (3-9) 


The  STY  objective  function  (Eq.  3-2)  is  calculated  as  the  NPV  of  a  terminating 
periodic  series  of  harvest  revenues  with  harvest  horizon  T=  40  years  and  cutting-cycle 
f  =  5  years,  which  is  the  most  profitable  interval  for  low-intensity,  sustained  harvesting 
given  the  effects  of  discounting.  Higher  NPV  is  achieved  at  shorter  cutting-cycles  for 
STY  and  SFM,  given  the  important  effects  of  an  1 8.75%  discount  rate.  Comparisons  of 
returns  from  harvesting  at  longer  cutting-cycles  were  conducted  in  preliminary  analyses 
to  confirm  the  validity  of  this  expectation.  The  sustainable  harvest  level  increases,  but 
the  present  value  of  net  returns  decreases  markedly  with  increases  in  cutting-cycle  length. 

Forest  growth  serves  as  the  sustained-yield  constraint  for  the  STY  case,  as  it 
requires  a  constant  periodic  harvest  (h,  =  hi+ys  =h*)  and  maintenance  of  a  specified  stand 
structure  (y,  =y,*-,o  =y*)-  Definition  of  the  sustained-yield  constraints  follows  that  of 
Buongiorno  and  Michie  (1980),  who  described  the  fundamental  mathematical  expression 
of  the  sustained-yield  linear  optimization  problem 

Grh'  +  {l-Gr)y='£Gkr  (3-10) 

*=o 

The  STY  regime  is  further  subject  to  constraints  defining  merchantable  proportion  (Eq. 
3-4)  and  size  (Eq.  3-6). 
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The  sustainable  forest  management  alternative  (SFM)  has  as  its  objective  to 
achieve  a  sustained-yield,  steady-state  forest  as  "close"  as  possible  in  structure  and 
composition  to  the  undisturbed  climax  condition.  The  objective  function  is  defined  as 
minimizing  the  sum  of  the  absolute  differences  between  the  number  of  trees  per  size  class 
and  guild  in  the  climax  forest  (ya)  and  in  the  SFM  residual  stand  (y**).  Mathematically, 
the  objective  function  is  expressed 

mm  V^.-v->-ZZ('VA+^)  (3-H) 

'    j 

where  [dj/]  and  [d,/]  are  vectors  of  the  absolute  deviations  of  the  number  of  trees  of  size 
i  and  guild/  in  the  growing  stocky**  and  the  climax  state/*  (Buongiorno  et  al.  1995) 
and  Bij  is  the  mean  basal  area  per  size  class  (/)  and  guild  (,),  employed  to  place  more 
weight  on  large  trees. 

The  SFM  alternative  observes  the  same  constraints  as  the  STY  case,  notably 
sustained-yield  forest  growth  (Eq.  3-10)  and  merchantability  restrictions  (Eq.  3-4,  3-6), 
and  is  further  constrained  to  meet  the  minimum  economic  goals  of  obtaining  a  non- 
negative  NPV  in  perpetual  periodic  harvest  (Eq.  3-12). 
v'ti'-F     A. 


(1  +  rf 


->0  (3-12) 


Residual  Stand  State  Indices 

Stand  state,  forest  value  and  diversity  criteria  are  employed  to  examine  changes  in 
the  managed  stand  under  alternative  regimes.  The  criteria  allow  for  a  comparison  of 
residual,  managed  stands  under  optimal  solutions  to  the  four  management  alternatives 
examined  and  comparison  of  these  states  relative  to  the  climax  condition  defined  by  the 
growth  model  (y").  Stand  basal  area  (BA)  is  a  fundamental  measure  of  change  in  stand 
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density  (m  /ha),  which  reveals  both  the  level  of  harvest  intensity  and  the  "openness"  of 
the  residual  stand  in  terms  of  soil  occupancy  by  trees.  Changes  in  structure  are  examined 
by  graphically  comparing  the  residual  distribution  of  trees  by  size  class  and  by  estimating 
the  proximity  of  residual  stand  structure  and  composition  to  the  climax  state.  The 
proximity  to  climax  index  (PCI)  introduced  by  Boscolo  and  Vincent  (2000)  is  employed 
to  examine  the  deviation  in  structure  and  composition  of  managed  stands  from  the  climax 
equilibrium.  Proximity  to  climax  structure  and  density  at  time  t  is  defined 


Et/ivM? 


";  zow  (3"I3) 

PCI  measures  the  proximity  of  the  residual  stand  population  at  time  /  (y,,,)  to  the  climax 
distribution  (y°  ,,).  The  mean  basal  area  per  size  class  and  guild  (By)  is  used  to  place  more 
weight  on  the  large  trees  that  are  associated  with  climax  forest  conditions  and  high 
endemic  biodiversity  (cf.,  Boscolo  and  Vincent  2000). 

In  addition  to  the  PCI  measure  of  deviation  from  climax  structure  and  density,  the 
PCI  form  is  used  to  examine  two  other  characteristics,  proximity  to  climax  market  value 
(PCM)  and  proximity  to  climax  wildlife  value  (PCW).  Structural  differences  are 
captured  in  PCI  and  thus  are  not  repeated  in  these  indices;  instead,  the  residual  density  of 
stems  per  sub-guild  (yhBk)  is  compared  to  the  climax  condition  (y"kBk).  Sub-guilds  (k) 
are  derived  from  Guilds  (J)  by  differentiation  of  merchantable  and  non-merchantable 
species,  except  for  Guild  3  that  contains  no  merchantable  species.  Proximity  to  climax 
market  value  at  time  t  (PCM,)  is  defined 


iZW'Wt-y^Mtf 


PCM.  -1-  M — ■* (3-14) 
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where  Mk  is  the  mean  ordinal  market  value  of  sub-guild  k  (Table  3-4),  derived  from  a 
ranking  of  tree  species  per  market  value  (0  =  no  value,  3  =  high),  which  is  a  refinement  of 
initial  rankings  by  Pinard  et  al.  (1999)  conducted  by  the  author  and  BOLFOR  specialists 
(Appendix  A).  The  PCM  measures  the  commercial  timber  stock  value  of  residual  forest 
states  (y,)  relative  to  that  of  the  climax  forest  (ya). 

Similarly,  proximity  to  climax  wildlife  value  (PCW)  is  defined 


PCW,=\-     '      _     a (3-15) 

where  Wk  is  the  mean  ordinal  wildlife  value  of  sub-guild  k  (Table  3-4),  derived  from  the 

ranking  of  tree  species  per  wildlife  value  (1  =  low,  3  =  high)  by  Pinard  et  al.  (1999)  and 

updated  with  BOLFOR  specialists  (Appendix  A).  Wildlife  values  are  ranked  by  the 

importance  of  tree  species  for  vertebrate  frugivores,  based  on  interviews  with  local 

hunters  and  stomach  contents  of  wild  game  (Guinart  1997,  Pinard  et  al.  1999).  PCW 

permits  the  interpretation  of  alternative  residual  stand  states  as  forest  habitat  for  resident 

wildlife. 

Table  3-4.  Mean  ordinal  market  and  wildlife  values  per  sub-guild 

Sub-guild Market  value  (Mk) Wildlife  value  (Wk) 

Guild  1,  merchantable  1.00  1.50 

Guild  1 ,  non-merchantable  0.00  2.16 

Guild  2,  merchantable  1.60  2.40 

Guild  2,  non-merchantable  0.00  2.33 

Guild  3,  non-merchantable  0.00  1.68 

Guild  4,  merchantable  2.33  1.00 

Guild  4,  non-merchantable  0.00  1 .57 

Guild  5,  merchantable  1.60  1.40 

Guild  5,  non-merchantable 0.00  1.77 
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The  Shannon- Wiener  index  (H),  with  natural  logarithms  (Whittaker  1975),  is  also 
used  to  examine  species  and  structural  diversity  as  a  measure  of  residual  stand  state, 
defined  as  follows: 

#»*»= -££>,>(#,)  (3-16) 

where  pr>  is  the  proportion  of  trees  in  the  size  class  i  and  guild/  For  the  Shannon-Wiener 
diversity  analysis,  trees  are  grouped  two  ways:  (1)  for  guild  diversity,  the  five  guild 
groupings,  and  (2)  for  guild  and  structural  diversity,  80  groups  defined  by  guild  and  size 
class  (dbh).  H  is  a  measure  of  evenness,  which  is  maximized  if  trees  are  equally 
distributed  by  size  class  and  guild,  //has  a  minimum  value  of  zero  for  a  homogeneous 
stand  and  maximum  value  of  ln(n+m)  for  a  perfectly  heterogeneous  stand;  consequently, 
the  ranges  of  possible  values  for  the  diversity  criteria  described  above  are  0  <  H5  <  1 .61 
for  guild  diversity  and  0  <  rf°<  4.38  for  guild  and  structural  diversity. 

Lastly,  the  absolute  difference  between  the  number  of  trees  per  size  class  and  guild 
in  the  climax  forest  (y°)  and  in  the  residual  stand  (y<)  defined  as  the  objective  function  of 
the  SFM  case  (Eq.  3-1 1)  is  used  to  examine  differences  in  alternative  managed  stands 
relative  to  the  climax  condition. 

Results  and  Discussion 
Optimal  Harvest  per  Management  Alternative 

Given  the  high  discount  rates  applicable  to  capita]  investments  in  Bolivia  and  the 
low  rate  of  value  growth  of  timber  species,  which  is  estimated  at  a  maximum  of  1 .23% 
assuming  constant  prices,  the  optima]  solution  for  an  unconstrained  harvest  scenario  (U) 
is  to  harvest  all  the  merchantable  timber  and  invest  the  returns  from  harvest  in  more 
lucrative  activities.  Consequently,  any  retention  of  merchantable  timber  would  result  in  a 
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decrease  in  the  profitability  of  management.  This  prescription  for  liquidation  of  the 

merchantable  timber  stock  is  found  to  be  the  optimal  solution  for  scenario  U,  in  which 

harvest  is  solely  constrained  by  the  merchantability  limits  established  in  local  timber 

markets.  Consequently,  harvest  level  is  highest  under  unconstrained  management  (Table 

3-5)  and  forest  impacts  are  greatest. 

Table  3-5.  Harvest  level  in  absolute  terms  and  in  proportion  of  the  available 
merchantable  volume  in  Year  0  per  management  alternative 


Case 

Harvest  level  (stems/ha) 
Y      .          YearO         Year  40, 
proportion    cumulative 

] 
YearO 

harvest  volume  (m3/hal 
Year  0         Year  40, 
proportion    cumulative 

U 

BFL 

STY 

SFM 

25.67 

13.36 

4.08 

0.51 

100% 
52% 
16% 

2% 

35.53 
19.91 
36.76 

4.57 

35.86 

24.71 

4.18 

1.63 

100% 
69% 
12% 

5% 

45.61 
33.75 
37.63 
14.67 

Harvest  constrained  by  the  Bolivian  forestry  law  has  a  similar  optimal  solution,  notably 
that  of  removing  all  merchantable  timber  allowed  by  legal  and  merchantable  limits.  The 
lower  harvest  level  and  inferior  net  returns  of  BFL  relative  to  the  unconstrained  case  are 
imposed  by  dbh  limits  and  the  allowable  cut  restriction,  which  requires  retention  of  20% 
of  the  merchantable  timber.  The  reduction  in  net  revenues  for  the  BFL  scenario 
attributable  to  the  45  cm  dbh  limit  and  80%  allowable  cut  constraints  on  the  initial 
harvest  entry  are  $145.45/ha  and  $191.60fta,  respectively.  The  optimal  harvest  solution 
for  BFL  prescribes  progressive  depletion  of  the  merchantable  stock,  subject  to  these  legal 
restrictions. 

The  optimal  sustained-yield  (STY)  solution  requires  a  substantial  decrease  in 
harvest  level  for  the  initial  entry  relative  to  U  and  BFL,  but  the  total  cumulative  volume 
harvested  over  a  40-year  management  horizon  is  greater  than  the  BFL  alternative  (Table 
3-5).  This  is  attributable  to  the  fact  that  STY  is  not  constrained  to  meet  dbh  limit  or 
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proportional  harvest  restrictions  as  is  BFL  and  to  the  expectation  that  the  STY  system 
better  conserves  the  merchantable  growing  stock  over  the  harvest  horizon.  The  STY 
management  for  a  5-year  cutting-cycle  stipulates  intensive  harvesting  of  stems  of 
merchantable  size  in  Guilds  4  and  5,  the  upland  forest  species  and  highly  selective 
harvesting  of  the  riparian  species  of  Guilds  1  and  2  (Table  3-6).  The  optimal  STY 
solution  is  determined  by  three  factors:  (1)  biological  constraints  of  forest  growth  and 
recruitment  and  model;  (2)  cutting-cycle  length;  and,  (3)  timber  merchantability. 
Changes  in  management  costs  and  discount  rates  have  no  effect  on  the  STY  solution. 


Table  3-6. 

Comparison  of  pre-harvest  population 
distribution  (y*)  per  guild 

(yo)  and  sustained 

-yield  (STY) 

Dbh  class 

Pre-harvest  population 

Vo  (trees/ha) 

STY  distribution,  v* 

(trees/hal 

(cm) 

Guild  1 

2 

3 

4 

5 

1 

2 

3 

4 

5 

10 

47.1 

2.7 

6.2 

23.0 

91.2 

38.8 

4.6 

6.6 

35.1 

109.5 

15 

28.3 

2.6 

5.0 

22.2 

83.1 

21.4 

3.4 

4.2 

29.3 

73.8 

20 

12.3 

l.X 

1.3 

17.2 

50.2 

12.9 

2.7 

2.7 

23.7 

49.7 

25 

5.8 

2.0 

0.5 

13.3 

26.2 

8.4 

2.2 

1.7 

18.2 

33.2 

30 

4.6 

0.9 

0.2 

10.9 

16.8 

5.8 

1.8 

1.1 

13.3 

21.9 

35 

3.1 

0.9 

0.4 

11.2 

10.8 

4.2 

1.5 

0.7 

9.1 

14.2 

41) 

2.8 

0.9 

0.5 

9.4 

6.2 

3.1 

1.3 

0.4 

2.1 

4.9 

45 

2.2 

0.8 

0.2 

5.4 

3.7 

2.3 

1.1 

0.2 

0.0 

1.2 

50 

1.6 

0.3 

0.3 

5.1 

1.9 

1.8 

1.0 

0.1 

0.0 

0.3 

55 

0.7 

0.4 

0.1 

3.2 

1.0 

1.4 

0.8 

0.1 

0.0 

0.1 

60 

1.0 

0.2 

0.1 

1.6 

0.7 

1.1 

0.7 

0.0 

0.0 

0.0 

65 

0.4 

0.2 

0.2 

0.6 

0.7 

0.8 

0.6 

0.0 

0.0 

0.0 

70 

0.3 

0.1 

0.2 

0.5 

0.2 

0.6 

0.5 

0.0 

0.0 

0.0 

75 

0.1 

0.0 

0.1 

0.2 

0.1 

0.5 

0.4 

0.0 

0.0 

0.0 

80 

0.2 

0.0 

0.1 

0.4 

0.0 

0.1 

0.3 

0.0 

0.0 

0.0 

85 

0.2 

0.1 

0.6 

0.1 

0.1 

0.0 

0.1 

0.0 

0.0 

0.0 

A  dramatic  change  in  the  forest  structure  is  prescribed  for  the  STY  condition 
relative  to  the  initial  stand  state,  especially  for  upland  species  of  Guilds  4  and  5  (Table  3- 
6).  With  frequent  harvest  entries,  STY  prescribes  maintaining  a  small  growing  stock  to 
generate  timber  just  sufficient  to  meet  sustained-yield  goals  and  generate  modest  profits. 
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Comparison  with  the  pre-harvest  stand  state  reveals  that  the  STY  condition  may  not  be 

assumed  from  Year  0,  rather  a  period  of  regulation  harvest  is  required  to  convert  the 

stand  to  the  STY  steady-state.  The  duration  of  the  regulation  period  and  the  returns  to 

management  aimed  at  achieving  the  STY  condition  are  unknown.  Estimating  the  returns 

to  regulation  harvest  is  not  as  simple  as  subtracting  the  STY  or  SFM  populations  from  the 

initial  population  and  estimating  the  value  of  the  implied  harvest.  Rather,  it  is  expected 

that  regulation  would  require  harvesting  both  merchantable  and  non-merchantable  stems 

and  active  management  toward  the  desired  STY  or  SFM  conditions  relative  to  growth 

and  regeneration  responses  to  silvicultural  interventions.  Thus,  while  it  may  be  expected 

that  the  returns  to  STY  or  SFM  would  be  increased  from  an  initial  harvest  intended  to 

convert  the  stand  to  desirable  conditions,  the  nature  and  value  of  these  returns  are 

uncertain. 

Table  3-7.  Comparison  of  pre-harvest  population  (yo)  and  sustained-yield  (SFM) 

distribution  (y**)  per  guild 

Dbh  class    Pre-harvest  population.  Vn  (trees/ha)         SFM  distribution,  v**  (trees/ha) 
(cm)     Guild  1      2  3  4  5  12  3  4  5 


10 

47.1 

2.7 

6.2 

23.0 

91.2 

38.8 

4.6 

6.6 

35.1 

109.5 

15 

28.3 

2.6 

5.0 

22.2 

83.1 

21.4 

3.4 

4.2 

29.3 

73.8 

20 

12.3 

1.8 

1.3 

17.2 

50.2 

12.9 

2.7 

2.7 

23.7 

49.7 

25 

5.8 

2.0 

0.5 

13.3 

26.2 

8.4 

2.2 

1.7 

18.2 

33.2 

30 

4.6 

0.9 

0.2 

10.9 

16.8 

5.8 

1.8 

1.1 

13.3 

21.9 

35 

3.1 

0.9 

0.4 

11.2 

10.8 

4.2 

1.5 

0.7 

9.1 

14.2 

40 

2.8 

0.9 

0.5 

9.4 

6.2 

3.1 

1.3 

0.4 

5.7 

9.0 

45 

2.2 

0.8 

0.2 

5.4 

3.7 

2.3 

1.1 

0.2 

3.4 

5.6 

50 

1.6 

0.3 

0.3 

5.1 

1.9 

1.8 

1.0 

0.1 

1.8 

3.3 

55 

0.7 

0.4 

0.1 

3.2 

1.0 

1.4 

0.8 

0.1 

0.8 

2.0 

60 

1.0 

0.2 

0.1 

1.6 

0.7 

1.1 

0.7 

0.0 

0.1 

1.1 

65 

0.4 

0.2 

0.2 

0.6 

0.7 

0.8 

0.6 

0.0 

0.0 

0.4 

70 

0.3 

0.1 

0.2 

0.5 

0.2 

0.6 

0.5 

0.0 

0.0 

0.1 

75 

0.1 

0.0 

0.1 

0.2 

0.1 

0.5 

0.4 

0.0 

0.0 

0.0 

80 

0.2 

0.0 

0.1 

0.4 

0.0 

0.4 

0.3 

0.0 

0.0 

0.0 

85 

0.2 

0.1 

0.6 

0.1 

0.1 

0.1 

0.1 

0.0 

0.0 

0.0 
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The  SFM  condition  prescribes  maintenance  of  the  negative  exponential  distribution 
typical  of  uneven-aged,  mixed  species  forests,  with  more  modest  reduction  in  the  pre- 
harvest  growing  stock  relative  to  STY  (Table  3-7).  The  SFM  condition  is  defined  by  the 
biological  needs  of  meeting  sustained  timber  yield  constraints  and  minimal  deviation 
from  the  theoretical  climax  structure.  Some  substitution  of  species  of  lower  market  value 
for  those  of  higher  value  may  occur  in  satisfying  the  minimum  deviation  objective,  due  to 
inadequacies  of  the  SFM  optimization  model.  Evidence  of  substitution  for  large  stems  of 
Guild  4,  those  of  highest  value,  is  inconclusive  given  the  important  stocking  levels 
observed  (Table  3-7). 

The  sustained-yield  scenarios  are  advantageous  in  that  the  harvest  horizons  are 
infinite,  thus  profitable  timber  extraction  may  be  concentrated  in  a  permanent  production 
forest  estate.  Optimal  unconstrained  management  (U)  depletes  the  merchantable  stock  in 
the  initial  entry  (Year  0).  The  BFL  alternative  only  retains  sufficient  merchantable  stock 
to  meet  legal  constraints,  implying  that  only  a  proportion  of  those  stems  growing  into  size 
classes  >  45  cm  during  a  cutting-cycle  may  be  harvested  in  subsequent  entries.  This 
result  implies  that  after  the  initial  concession  period  (40  years),  U  and  BFL  would  require 
a  growing  interval  longer  than  40  years  for  the  forest  to  produce  adequate  timber  for 
profitable  harvesting.  Removing  the  stand  from  production  for  several  decades  to  allow 
for  sufficient  restocking  of  merchantable  timbers  after  U  and  BFL  harvesting  creates  a 
short-term  need  for  the  expansion  of  harvesting  into  other  forests  to  meet  production 
goals,  which  defeats  the  purpose  of  establishing  permanent  production  forests. 
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Profitability  of  Management 

Economic  returns  from  the  alternative  management  scenarios  indicate,  as  expected, 
a  considerable  decrease  in  the  NPV  of  sustained-yield  harvesting  relative  to 
unconstrained  harvesting  (Table  3-8).  The  opportunity  costs  of  the  sustained-yield 
alternatives  (STY,  SFM)  are  likewise  very  high,  exceeding  $1000  per  hectare  for  both 
scenarios  relative  to  unconstrained  liquidation  of  the  available  merchantable  timber  stock, 
assuming  accuracy  of  the  cost  assumptions.  The  comparative  analysis  of  financial 
returns  furthermore  assumes  that  there  is  no  period  of  regulation  harvesting,  thus  no  costs 
or  financial  benefits  derived  from  converting  the  forest  from  its  initial  state  (yo)  to  the 
STY  (y*)  or  SFM  (y**)  conditions. 

The  great  disparity  in  returns  from  sustained-yield  and  unconstrained  alternatives  is 
attributable  to  the  opportunity  costs  of  leaving  merchantable  stems  in  initial  harvest 
entries,  given  the  high  discount  rate  (18.75%).  High  opportunity  costs  of  sustained-yield 
management  would  occur  at  even  lower  discount  rates,  however,  given  the  low  value 
growth  rate  of  the  tropical  hardwoods  in  this  Chiquitano  forest. 

Table  3-8.  Net  present  value  and  opportunity  costs  of  management  alternatives 

NPV  Opportunity 

Case ($/ha) costs  ($/ha) 

U  1380.54  0.00 

BFL  914.89  465.64 

STY  257.18  1123.36 

SFM 1L25 1369.29 

The  sustained-yield  scenario  (STY)  yields  profits  of  $257.18/ha  despite  the  high 
discount  rate.  The  more  conservative  SFM  harvest  scenario  also  yields  profits 
($1 1 .25/ha)  while  more  closely  approaching  the  climax  condition.  Optimal  management 
under  constraints  imposed  by  Bolivia's  Forestry  Law  (BFL)  yields  superior  financial 
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returns  relative  to  the  sustained-yield  alternatives,  but  has  important  opportunity  costs 

relative  to  unconstrained  harvesting  ($465.64/ha).  These  opportunity  costs  derive  from 

merchantable  stems  not  harvested  in  Year  0  in  order  to  meet  dbh  and  allowable-cut 

limits. 

The  discount  rate  affects  the  profitability  of  management  under  those  alternatives 

whose  objective  function  is  to  maximize  this  output  (U,  BFL,  STY),  but  does  not  change 

the  optimal  solution.  The  optimal  solution  is  constrained  by  biological  and 

merchantability  restrictions  and  cutting-cycle  length.  At  lower  discount  rates,  STY 

becomes  more  competitive  with  the  U  and  BFL  options  (Table  3-9),  exceeding  the 

profitability  of  BFL  under  1%  discounting.  The  SFM  optimal  harvest  solution  changes 

with  the  discount  rate  as  less  deviation  from  the  climax  condition  is  allowed  at  lower 

discount  rates,  while  satisfying  the  binding  constraint  of  non-negative  NPV  for  a 

perpetual  harvest  horizon.  Consequently,  the  NPV  of  SFM  harvesting  for  a  40-year 

horizon  decreases  with  decreases  in  the  discount  rate. 

Table  3-9.  Discount  rate  effects  on  the  NPV  of  harvest  returns  for  a  40-year  horizon 

NPV  ($/ha) 


Case\  Discount  rate: 

18.75% 

12.00% 

6.00% 

1.00% 

U 

1380.54 

1379.72 

1402.11 

1593.29 

BFL 

914.89 

913.55 

931.31 

1091.75 

STY 

257.18 

338.98 

538.41 

1087.67 

SFM 

11.25 

9.84 

8.73 

7.90 

Stand  State  and  Forest  Value  Measures 

Comparison  of  the  structure  of  aggregate  merchantable  and  non-merchantable 
residual  populations  under  the  alternative  management  regimes  and  that  projected  for  the 
undisturbed  climax  reveals  minor  differences  in  stand  impact  among  management  options 
(Fig.  3-1).  As  expected,  the  unconstrained  scenario  has  the  greatest  impact  on  residual 
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stand  structure,  as  the  merchantable  timber  stock  is  progressively  liquidated.  Also  as 
anticipated,  the  SFM  regime  closely  approximates  the  climax  structure,  meeting  its 
defined  objective  of  minimizing  deviation  from  this  undisturbed,  equilibrium  state. 


U,  40  yrs 

BFL,  40  yrs 

STY  condition 

SFM  condition 

Climax 


25   30   35   40   45   50   55   60   65   70   75   80   85 

Dbh  class  (cm) 

Figure  3-1.  Aggregate  merchantable  and  non-merchantable  residual  population  structure 
under  management  alternatives  and  climax  condition 

Aggregate  residual  population  structures  do  not  differ  greatly  under  U  and  STY 

management  regimes,  however,  and  the  BFL  constraints  appear  to  achieve  greater 

proximity  to  climax  stand  structure  than  the  STY  alternative.  Proximity  to  climax 

structure  and  composition  of  the  aggregate  residual  populations  (PCI)  confirms  that  the 

BFL  alternative  results  in  greater  proximity  to  this  undisturbed  forest  equilibrium  than  U 

and  STY,  as  indicated  in  the  comparison  of  aggregate  population  structure  (Table  3-10). 

This  finding  is  unexpected,  but  may  be  explained  by  the  fact  that  the  aggregate 

population  may  be  effectively  maintained  under  BFL,  while  the  stand  is  impoverished. 
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This  finding  is  unexpected,  but  may  be  explained  by  the  fact  that  the  aggregate 
population  may  be  effectively  maintained  under  BFL,  while  the  stand  is  impoverished. 


Table  3-10. 

State  of  the 

aggregate  residual  stand  populatior 

i  per  management  alternative 

Management 

alternative 

Climax 

Index 

u 

BFL 

STY 

SFM 

condition 

BA 

19.72 

20.76 

21.72 

26.19 

26.90 

PCI 

0.73 

0.78 

0.68 

0.97 

1.00 

D 

7.50 

6.46 

5.18 

0.56 

0.00 

This  hypothesis  is  confirmed  in  closer  examination  of  the  impacts  of  harvest  on  the 

merchantable  stock  (Table  3-1 1). 

Table  3-11.  Residual  merchantable  species  population  per  management  alternative 

Management  alternative Climax 


Index 

U 

BFL 

STY 

SFM 

condition 

BA  (all  merch.  species) 

5.42 

6.46 

8.55 

9.83 

12.48 

BA  (merch.  species  >  40  cm  dbh) 

0.00 

1.04 

2.64 

3.92 

5.57 

"^tmcrch 

0.42 

0.53 

0.57 

0.76 

1.00 

PCM 

0.62 

0.70 

0.83 

0.86 

1.00 

Note:  PCImerch  is  defined  as  PCI,  but  calculated  solely  for  merchantable  species. 

Residual  stand  basal  area  (BA)  measures  reveal  greater  differences  between  the  U 
and  STY  alternatives  relative  to  comparisons  of  the  aggregate  residual  populations, 
though  the  differences  of  these  states  and  BFL  are  less  important  (Table  3-10).  Under 
timber  management  objectives  of  U,  BFL  and  STY,  residual  stand  density  is  expected  to 
be  19%  to  27%  lower  than  the  climax  forest.  Harvesting  impacts  on  the  merchantable 
stock  are  more  revealing  of  differences  between  the  management  alternatives,  as 
suggested  by  the  proximity  to  climax  market  value  (PCM)  index  (Table  3-1 1).  The  PCM 
of  residual  stands  is  considerably  lower  under  U  and  BFL  regimes,  than  under  STY, 
given  the  effective  liquidation  of  merchantable  stems  stipulated  by  these  scenarios 
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prescribes  maintenance  of  a  growing  stock  of  merchantable  stems  for  future  harvest  (Fig. 
3-2). 

Impacts  of  the  alternative  harvest  regimes  on  the  merchantable  species  stock  are 
further  elucidated  in  examination  of  residual  basal  area  (BA),  proximity  to  climax 
merchantable  stocking  (PCImerch).  and  proximity  to  climax  market  value  (PCM)  of  the 
residual  stands  (Table  3-1 1).  These  indices  reveal  much  greater  impacts  of  U  and  BFL 
regimes  over  the  concession  agreement  (40  years).  The  basal  area  comprised  by 
merchantable  species  of  all  size  classes  under  unconstrained  harvesting  is  reduced  to  less 
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Figure  3-2.  Residual  merchantable  population  structure  under  management  alternatives 
and  climax  condition 

than  half  that  of  the  climax  forest  and  to  just  over  half  of  the  climax  merchantable 

stocking  under  BFL  constraints. 
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than  half  that  of  the  climax  forest  and  to  just  over  half  of  the  climax  merchantable 
stocking  under  BFL  constraints. 

The  basal  area  of  merchantable  stems  is  depleted  under  U  and  reduced  to  less  than 
1/5  of  that  in  the  climax  forest  under  BFL.  In  contrast,  the  STY  and  SFM  regimes 
maintain  important  stocking  levels  of  merchantable  stems  throughout  the  harvest  horizon. 
Moreover,  significant  differences  are  observed  in  the  proximity  to  climax  structure  and 
composition  indices  (PCImerch)  under  scenarios  constrained  to  sustained  timber  yield 
(STY,  SFM)  and  those  not  biologically  constrained  to  meet  sustainability  objectives 
(U,  BFL).  The  PCM  measures  confirm  the  depletion  of  stand  market  value  under  the  U 
and  BFL  relative  to  the  sustained-yield  alternatives. 

Wildlife  values  of  the  residual  stands  are  less  affected  by  the  alternative 
management  scenarios  (Table  3-12),  indeed  all  scenarios  result  in  residual  stands  that 
deviate  little  from  the  climax  stand  in  their  value  to  vertebrate  frugivores.  This  result  is 
expected,  as  non-merchantable  species  are  as  important  as  merchantable  species  for 
wildlife,  if  not  more  so  (Table  3-4).  Accordingly,  depletion  of  the  merchantable  stock 
has  little  impact  on  the  value  of  the  forest  for  wildlife,  as  measured  in  this  study.  Greater 
precision  on  the  relative  value  of  tree  species  to  wildlife  other  than  vertebrate  frugivores 
may  be  more  revealing  of  the  impacts  of  alternative  management  scenarios  in  terms  of 
forest  value  for  biodiversity. 
Table  3-12.  Wildlife  and  diversity  indices  per  management  alternative 


Management 

alternative 

Climax 

Index 

U 

BFL 

STY 

SFM 

condition 

PCW 

0.98 

0.99 

0.96 

0.98 

1.00 

H5 

1.17 

1.17 

1.21 

1.20 

1.20 

Hg0 

2.93 

2.97 

2.98 

3.04 

3.11 
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Shannon- Wiener  measures  of  stand  diversity  also  reveal  little  difference  among 
management  alternatives  (Table  3-12).  The  SFM  results  in  slightly  greater  evenness 
across  guilds  and  size  classes  than  would  be  expected  in  under  more  intensive  harvest 
scenarios  (U,  BFL,  STY),  as  harvesting  results  in  a  decrease  in  merchantable  stems.  The 
Shannon- Wiener  diversity  indices  reveal  little  about  differences  among  these  timber- 
focused  management  alternatives. 

Conclusions 

Management  of  the  Chiquitano  dry  tropical  forests  of  eastern  Bolivia  for  sustained- 
timber  yield  is  shown  to  be  feasible  and  profitable,  given  certain  simplifying  assumptions 
of  the  forest  growth  model.  Moreover,  the  present  study  demonstrates  that  sustained- 
yield  management  with  greater  conservation  of  forest  structure  and  composition  may  be 
profitably  implemented.  The  optimal  trade-off  between  STY  and  SFM  conditions  is  not 
clearly  defined  and  will  depend  upon  further  information  concerning  the  marginal 
ecological  impacts  of  deviation  from  the  theoretical  climax  forest  condition.  The  optimal 
sustained  timber  yield  (STY)  condition,  defined  as  one  that  maximizes  the  returns  to 
harvest  given  sustained  production  constraints,  is  determined  by  biological  criteria  for 
sustainability,  cutting-cycle  length  and  timber  merchantability  constraints.  Harvesting 
costs  and  discount  rates  have  no  effect  on  the  optimal  STY  solution,  unless  costs  increase 
to  the  point  of  rendering  harvest  unprofitable,  which  would  eliminate  consideration  of 
logging  activities.  Marginal  changes  in  costs  and  discount  rates  have  not  effect  on  the 
optimum  given  the  important  constraints  posed  by  the  biological  requirements  for 
sustainability. 

Forest  management  under  STY  and  SFM  steady-state  conditions  is  assumed 
possible  in  Year  0  to  allow  for  comparative  analysis  of  these  theoretical  forest 
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management  ideals  with  conventional  harvest  options  (U  and  BFL).  This  presumes  that 
the  regulation  harvest  necessary  to  convert  the  forest  to  sustained-yield  conditions  is 
feasible.  Moreover,  the  assumption  of  steady-state  conditions  implies  that  the  duration 
and  magnitude  of  regulation  harvest  do  not  influence  the  returns  to  management.  This 
assumption  is  untenable  for  management  of  the  Las  Trancas  forests  from  their  current 
state;  however,  strict  comparison  of  returns  to  the  theoretical  STY  and  SFM  regimes  with 
conventional  alternatives  remains  useful  and  informative  of  the  financial  and  ecological 
implications  of  management  under  sustainable  equilibria.  While  it  may  be  expected  that 
the  returns  to  STY  or  SFM  would  be  increased  from  an  initial  harvest  intended  to  convert 
the  stand  to  desirable  conditions,  the  nature  and  value  of  these  returns  are  uncertain. 
Estimating  the  management  requirements  and  duration  of  regulation  harvest  will  be  an 
important  extension  of  the  present  work. 

The  comparison  of  optimal  harvesting  behavior  with  and  without  constraints  aimed 
at  maintaining  the  productivity,  structure  and  composition  of  Chiquitano  dry  forests 
indicates,  as  anticipated,  that  the  opportunity  costs  of  more  sustainable  management 
relative  to  unconstrained  harvesting  are  considerably  high,  assuming  modest  increases  in 
profitability  from  conversion  to  the  STY  and  SFM  conditions.  The  magnitude  of  the 
decline  in  harvest  profitability  of  sustained  timber  yield  (STY)  relative  to  unconstrained 
(U)  management  may  be  largely  attributed  to  the  effects  of  excessively  high  discount 
rates  applicable  to  capital  investments  in  Bolivia.  Nonetheless,  given  the  low  value 
growth  rates  of  tropical  hardwoods  harvested  in  the  Chiquitano  forests,  the  requirement 
for  leaving  merchantable  timber  in  the  forest  imposed  by  both  the  STY  and  SFM  regimes 
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implies  that  these  options  will  consistently  generate  inferior  returns  from  unconstrained 
liquidation  of  merchantable  timber  in  initial  harvest  entries. 

Examination  of  the  economic  and  biological  impacts  of  an  optimal  management 
regime  constrained  by  the  Bolivian  forestry  law  indicates  that  although  financial  returns 
from  management  under  BFL  are  robust,  harvesting  constraints  are  not  likely  to  achieve 
goals  of  sustainable  timber  production  or  ecosystem  maintenance.  This  conclusion 
cannot  be  drawn  without  caveats,  given  the  simplifying  assumptions  required  for 
modeling  the  BFL  scenario  with  the  available  growth  and  yield.  Greater  precision  in 
modeling  diameter  limits  and  allowable  cut  restrictions  for  individual  species  and  in 
projecting  forest  growth  and  regeneration  may  indicate  that  optimal  BFL  management 
results  in  forest  conditions  more  closely  approximating  a  sustained-yield  state  than  is 
possible  to  project  in  this  study.  However,  the  simplified  model  used  in  the  present 
analysis  indicates  that  tree  diameter  and  allowable  cut  restrictions  lead  to  the  depletion  of 
large  stems  of  high  value,  which  is  not  likely  to  result  in  the  maintenance  of  forest 
structure  and  composition  necessary  to  achieve  sustained  timber  yield. 

Results  indicate  that  management  regimes  are  possible  that  achieve  sustained 
timber  yield,  while  generating  profitable  returns.  Moreover,  it  is  shown  that  maintaining 
forest  structure  and  composition  can  be  compatible  with  producing  sustained  timber  yield 
to  meet  economic  objectives  in  management  of  the  Chiquitano  dry  forests,  as  shown  by 
the  SFM  scenario.  For  the  particular  case  of  Chiquitania,  goals  of  sustainable  timber 
production  appear  compatible  with  the  maintenance  of  forest  wildlife  values,  though 
more  information  on  the  values  of  timber  species  for  other  taxa  would  allow  for  a  more 
conclusive  assertion  of  this  apparent  compatibility. 


CHAPTER  4 

REGULATORY  POLICY  EFFICIENCY  IN  PROMOTING  SUSTAINABLE  TIMBER 

MANAGEMENT  FOR  A  BOLIVIAN  TROPICAL  DRY  FOREST 

Introduction 

Regulatory  policies  are  intended  to  provide  a  corrective  influence  on  the  private  use 
of  forestland  to  properly  defend  public  goods,  both  in  maintaining  desirable  benefit  flows 
from  natural  resources  and  in  mitigating  negative  externalities.  Fiscal  regulatory 
instruments,  such  as  royalty  systems,  have  also  been  identified  as  a  critical  determinant  of 
the  manner  and  magnitude  of  tropical  forest  degradation  by  providing  incentives  for 
logging  practices  such  as  high-grading  and  conversion  to  alternative  land  uses  (Repetto 
and  Gillis  1988,  Vincent  1990,  Hyde  and  Sedjo  1992,  Van  Kooten  and  Bulte  2000). 

Despite  the  important  influence  of  government  policies  on  tropical  forest  use,  few 
rigorous  comparisons  exist  of  the  impact  of  alternative  regulatory  mechanisms  on  private 
revenues,  government  rent  capture,  and  forest  condition.  Important  exceptions  include 
recent  studies  by  Boscolo  and  Vincent  (2000)  and  by  Amacher  et  al.  (2001)  for  natural 
tropical  forests  of  Malaysia.  Boscolo  and  Vincent  (2000)  examined  concession  length, 
renewability  and  performance  bonds  as  mechanisms  to  promote  the  adoption  of  better 
logging  practices,  while  Amacher  et  al.  (2001)  used  a  model  of  policy  choice  to  compare 
royalty  systems  relative  to  government  revenue  generation  and  high-grading  behavior. 
Although  forest  degradation  has  often  been  attributed  to  poor  regulatory  policy,  there  is 
little  empirical  work  examining  how  selective  harvesting  behavior  and  the  allocation  of 
rent  between  private  concessionaires  and  government  is  influenced  by  regulatory 
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structure,  Amacher  et  al.  (2001)  being  the  exception.  The  present  study  aims  to 
contribute  to  such  empirical  analyses  by  considering  the  impacts  of  existing  regulatory 
policy  and  alternative  incentive  mechanisms  upon  optimal  harvesting  behavior  for  a 
tropical  dry  forest  of  eastern  Bolivia.  The  study  expands  upon  this  recent  work  by 
examining  the  influence  of  fiscal  regulatory  mechanisms  on  the  logging  behavior  of 
concessionaires,  the  private  returns  to  logging  and  government  rent  capture.  A  harvest 
optimization  model  is  used  to  conduct  a  comparative  analysis  of  the  efficiency  of 
alternative  policies  in  fostering  forest  management  compatible  with  sustained  timber 
yield  (STY)  objectives. 

Methods 
Forest  Growth  and  Optimization  Models 

The  seasonally  dry  tropical  forests  of  the  Lomerio  region  (16°13'S,  61°50'W)  lie  in 
a  transition  zone  between  the  humid  forests  of  the  southern  rim  of  the  Amazon  basin  and 
the  thorn  scrub  of  the  Gran  Chaco  on  the  southwestern  edge  of  the  Brazilian  shield 
(Killeen  et  al.  1998).  The  forests  are  typical  of  the  Chiquitano  dry  forests  of  eastern 
Bolivia.  Forest  growth  and  yield  is  projected  with  a  five-guild  matrix  growth  model 
estimated  by  multinomial  logistic  (MNL)  regression  of  permanent  sample  plot  data 
collected  in  two  400  ha  forest  blocks  near  the  Las  Trancas  community  (Claros  and  Licona 
1995,  Killeen  et  al.  1998).  The  MNL  matrix  growth  model  projects  the  probabilities  of 
upgrowth,  stability,  or  mortality  of  trees  in  each  of  five  ecological  guilds  during  a  5-year 
growth  period  and  predicts  recruitment  of  trees  into  the  smallest  size  class  (10  cm  dbh) 
for  each  guild  (cf.,  Chapter  2).  The  linear  form  of  the  MNL  model  and  the  assumption  of 
stationarity  of  transition  probabilities  throughout  the  harvest  horizon  enable  its 
straightforward  integration  into  this  linear  optimization  study. 
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The  Bolivian  forestry  law  was  revised  in  1996  (Ley  1700)  to  more  effectively 
promote  sustainable  management  of  the  nation's  forest  resources,  which  are  allocated  in 
renewable  concessions  to  the  logging  industry  (Art.  #1).  Logging  concessions  are 
awarded  for  a  harvest  period  of  40  years  upon  formal  application  to  and  approval  by 
Bolivia's  forest  service  (Superintendencia  ForestaT).  The  forestry  law  requires  forest 
management  plans  for  all  concessions  and  for  forest  authorizations  in  private  lands. 
Concessions  also  can  be  traded  and  inherited  and  are  renewable  upon  validation  of  the 
concessionaire's  observance  of  sustainable  forest  management  plans. 

Forest  management  costs  for  industrial  firms  of  the  region  and  market  prices  for 
roundwood  timber  delivered  to  the  forest  mill  derive  from  surveys  conducted  in  Santa 
Cruz  and  in  Chiquitania  in  2000  and  2001  (Appendix  B).  Prices  are  reported  in  2001 
$US  and  are  assumed  constant  throughout  the  harvest  horizon  (Appendix  C). 
Management  costs  in  2001  $US  are  classified  as  one  of  three  types: 

•  Variable  costs  (C)  in  $/m3  for  felling,  skidding  and  log  deck  operations  incurred 
relative  to  harvest  intensity 

•  Fixed  costs  (F)  in  $/ha  incurred  regardless  of  harvest  intensity  at  each  cutting-cycle 
entry  for  planning  and  capital  costs 

•  Annual  costs  (A)  in  $/ha  paid  throughout  the  harvesting  horizon  regardless  of 
harvest  intensity,  which  include  an  area  fee  of  $l/ha  for  the  concession  (Bolivia's 
patente)  and  inscription  fees  to  government  and  market  institutions 

These  and  all  other  variables  used  in  the  present  study  are  defined  in  Appendix  D. 

The  optimization  model  prescribes  maximization  of  net  present  value  (NPV)  of 

polycyclic  harvesting  (Z)  subject  to  biological,  market,  and  regulatory  constraints.  The 

objective  function  does  not  consider  the  costs  of  investment  in  the  growing  stock 

described  by  the  economic  stocking  rule  for  selection  harvest  (Duerr  and  Bond  1952)  or 

soil  expectation  value  (SEV)  used  in  similar  optimization  studies  (Buongiorno  and 
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Michie  1980,  Boscolo  and  Buongiorno  1997,  Bach  1999,  Kant  1999).  These  costs  are 
opportunity  costs  of  not  harvesting  merchantable  stems  in  the  initial  harvest  entry  in 
order  to  maintain  the  timber  growing  stock.  The  decision  to  exclude  the  opportunity 
costs  of  maintaining  the  growing  stock  from  the  calculation  of  net  returns  is  based  upon 
the  nature  of  the  forest  property.  As  these  forests  are  allocated  in  concessions  to  private 
industry  under  government  ownership,  the  opportunity  costs  of  conserving  the  residual 
stock  are  public  rather  than  private  costs. 

The  objective  function  (Z(h,  y>)  of  the  forest  management  problem  is  to  maximize 
the  NPV  of  returns  to  harvesting  over  a  40-year  concession  period  calculated  as  the 
present  value  of  a  terminating  series  at  discount  rate  8,  with  harvest  horizon  T=  40  years 
and  cutting  cycle  length  t  =  yd,  with  period  yd  being  the  number  of  intervals  (y)  of  5-year 
growth  periods  (0)  in  a  cutting-cycle.  For  the  present  analysis,  the  cutting  cycle  (t)  is  40 
years  and  y  =  8.  The  cutting-cycle  is  set  equal  to  the  length  of  the  concession  period 
because  the  returns  from  unconstrained,  profit  maximizing  harvesting  (U)  are  negative 
for  shorter  cutting-cycles  and  the  unconstrained  optimization  is  used  as  the  basis  for 
modeling  the  impact  of  regulatory  policies  (Chapter  3).  It  is  assumed  that  harvest  is 
allowed  in  Year  40,  before  the  end  of  the  concession  agreement.  The  objective  function 
is  defined  mathematically 

The  net  value  of  an  extracted  stem  of  size  class  i  in  Guild;  in  2001  $US/m3  (Table 
3-3)  is  calculated  as  the  product  of  the  mean  price  per  guild  (pj)  less  the  variable  costs  of 
harvesting  (C)  and  the  mean  extracted  volume  per  size  class  and  guild  (^,y).  Harvesting 
from  each  guild  (hj,)  is  measured  in  stems  per  ha.  The  tree  population  (y,)  and  harvest 
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intensities  (A,)  are  not  steady-state  variables,  as  the  regulatory  solutions  are  not  assumed 
to  achieve  sustained-yield  conditions.  The  real  discount  rate  (<5)  of  18.75  %  is  the  mean 
real  interest  rate  for  2001  reported  by  Bolivia's  Central  Bank  (Banco  Central  de  Bolivia 
2002). 

The  optimization  is  constrained  by  the  forest  growth  model  (Eq.  4-2),  which 
defines  growth  and  recruitment  in  the  residual  forest  between  harvest  intervals.  The 
forest  growth  constraint  requires  that  the  sum  of  the  harvest  and  growing  stock  at  time  t 
equal  the  residual  stocking  at  time  M  plus  any  growth  and  recruitment  in  this  growing 
stock  between  these  cutting-cycles. 

h,+y,-G'-(y,_l-h,J  =  fiGtr  (4-2) 

The  profit  maximization  objective  is  further  constrained  by  merchantability  criteria 

for  harvested  stems.  Merchantability  restrictions  require  that  the  proportion  of  trees 

harvested  per  Size  Class  ;'  and  Guild/  at  Time  t  be  equal  or  less  than  the  proportion  of 

merchantable  stock  (»/)  present  in  the  stand.  The  merchantability  constraint  for  the 

initial  harvest  (Year  0)  is  defined 

*„-«,.)>„  £0  (4-3) 

where  co,  is  the  merchantable  proportion  per  guild.  For  all  ensuing  cutting-cycles  under 

U  and  BFL,  the  merchantable  proportion  is  expected  to  decrease  as  merchantable  volume 

is  removed.  Thus  while  equation  4-3  is  applied  as  a  constraint  to  the  initial  harvest, 

subsequent  merchantable  stock  constraints  are  defined  by  the  growing  stock  of  residual 

merchantable  stems  remaining  after  harvest  (ft-i)  and  any  merchantable  upgrowth  into 

this  stock  during  the  cutting-cycle  interval.  This  constraint  is  described  mathematically 
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^.-G'jV.-IXgV^O  (4-4) 

In  addition,  merchantability  standards  require  that  trees  equal  or  exceed  40  cm  dbh: 

f>„,=0  (4-5) 

;-0 

Lastly,  the  optimization  is  subject  to  absolute  harvest  restrictions,  requiring  that  the 
harvest  intensity  not  exceed  the  growing  stock  (yt  -  h,  >  0)  and  non-negativity  constraints 
(ybht>0). 

Regulatory  Policies 

The  six  regulatory  policies  examined  are  as  follows: 

Bolivian  forestry  law  (BFL) 
Area  fee 
Per-tree  royalty 
Volume-based  royalty 
Ad  valorem  royalty 
Performance  bond 

The  level  of  fee  levied  for  each  policy  other  than  BFL  was  determined  iteratively, 

based  upon  the  target  of  achieving  a  residual  merchantable  population  approaching  that 

defined  for  the  STY  condition  (207  trees/ha  of  all  size  classes  and  guilds).  This  criterion 

was  selected  as  a  mechanism  for  examining  behavioral  change  relative  to  the  selective 

harvest  of  timber  species,  with  the  objective  of  comparing  the  efficiency  of  alternative 

regulations  in  fostering  the  maintenance  of  residual  forest  conditions  compatible  with  the 

structure  and  composition  necessary  for  sustained-yield  harvesting  (STY).  It  is  assumed 

that  the  administrative  costs  of  alternative  regulatory  mechanisms  are  equal.  Relaxing 

this  assumption  is  expected  to  affect  government  revenues  but  not  the  concessionaire's 

optimal  harvest  behavior  or  profits  (cf.,  Amacher  et  al.  2001).  Mathematical  descriptions 
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of  the  optimization  models  for  the  six  regulatory  policies  and  for  the  STY  model  are 
provided  in  Appendix  F. 

A  more  exhaustive  analysis  of  alternative  fee  levels  per  regulatory  mechanism 
would  enable  determination  of  the  optimal  fee  schedule  appropriate  to  achieve  both 
profitability,  rent  capture,  and  sustainability  goals.  This  analytical  approach  is  not  fully 
developed  in  the  present  study.  Instead,  fee  levels  required  to  compel  the  logger  to 
maintain  a  residual  merchantable  population  most  closely  approaching  that  required  for 
STY  were  selected  by  iterative  estimation  in  preliminary  model  simulations.  The 
alternative  regulatory  mechanisms  are  then  examined  at  the  selected  fee  levels  to  more 
closely  examine  behavioral  impacts  and  effects  on  financial  returns  from  logging  of  the 
policies,  without  detailed  examination  of  the  marginal  effect  of  alternative  fee  levels  per 
policy.  A  more  rigorous  examination  of  the  optimal  fee  levels  per  regulatory  mechanism 
would  be  a  valuable  extension  on  this  work,  though  the  present  analysis  provides  a  useful 
comparison  of  policy  impacts. 

Sustainable  timber  management  (STY).  The  regulatory  goal  of  this  analysis  is 
the  maintenance  of  forest  conditions  necessary  for  ST  Y;  to  wit,  maintaining  the 
productivity  of  the  forest  such  that  a  constant  harvest  yield  is  extracted  at  each  cutting- 
cycle  entry.  The  conditions  necessary  for  STY  are  derived  from  a  static  optimization 
analysis  using  the  MNL  growth  model  of  the  Las  Trancas  forest  for  a  cutting-cycle  of  5 
years  (cf.,  Chapter  3).  The  STY  harvest  intensity  is  estimated  at  approximately  4.2  m3/ha 
per  5-year  cutting-cycle. 

An  important  caveat  is  necessary  regarding  the  definition  of  the  forest  condition 
necessary  for  STY.  The  STY  condition  is  derived  from  optimization  analysis  using  the 


64 

MNL  forest  growth  model  described  in  Chapter  2.  The  optimal  STY  condition  is  largely 
dependent  upon  the  biological  constraints  imposed  by  this  growth  model,  by 
merchantability  constraints  for  timber  species  of  this  Chiquitano  forest  site,  and  by  the 
cutting-cycle  length  used  in  the  STY  analysis  (Chapter  3).  Alternate  STY  conditions  are 
possible  with  changes  in  these  fundamental  constraints.  The  present  study  assumes  that 
the  STY  condition  estimated  in  Chapter  3  is  a  desired  forest  condition  and  is  a 
sufficiently  accurate  definition  of  the  conditions  required  for  sustainability.  This  latter 
assumption  may  be  violated  for  long-term  modeling  of  sustainability  due  to  the  poor  fit  of 
the  model  of  forest  recruitment  (Chapter  2),  which  determines  long-term  changes  in 
forest  composition.  Therefore,  the  precise  conditions  necessary  for  STY  may  not  be 
conclusively  asserted  in  the  present  study.  Despite  this  important  weakness,  the 
comparative  analysis  of  the  effect  of  alternative  policy  mechanisms  on  logger  behavior 
remains  valid.  Its  validity  derives  from  the  fact  that  the  influence  of  regulatory 
mechanisms  is  not  dependent  upon  precise  definition  of  the  STY  condition. 

Bolivian  forestry  law  (BFL).  The  Bolivian  forestry  law  (Ley  1 700)  and  the 
technical  standards  for  forest  management  in  concessions  and  private  forests  greater  than 
200  ha  (Normas  Tecnicas)  established  normative  prescriptions  intended  to  promote  the 
maintenance  of  sustained  economic,  ecological,  and  social  benefits  of  production  forests. 
Under  this  legislation,  harvesting  is  restricted  to  the  following  standards:  a  minimum  20- 
year  cutting-cycle;  a  maximum  80%  allowable  cut  of  the  merchantable  stock  per  species; 
and,  a  minimum  dbh  limit  for  harvested  stems.  These  regulatory  constraints  are  imposed 
upon  the  optimization  model  to  examine  the  profitability  and  impact  of  expected 
management  behavior  under  the  Bolivian  forestry  law  in  the  "BFL"  scenario.  Under  the 
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BFL  case,  a  minimum  diameter  of  45  cm  is  required  for  merchantable  stems  (Eq.  4-6) 
and  an  80%  allowable  cut  restriction  is  applied  per  guild  (Eq.  4-7). 
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An  annual  area  fee  (palenle)  of  $l/ha  is  paid  on  all  logged  concession  lands  under 
terms  of  the  forest  law.  Concessionaires  may  set  aside  up  to  30%  of  the  forest  concession 
from  logging  and  be  exempted  from  paying  \hspatente  on  these  lands.  This  fee  is 
represented  in  the  objective  function  of  the  BFL  optimization  by  a  $l/ha  increase  in  the 
annual  costs  of  management  (A). 

Area  fee.  An  area  fee  of  (r  $/ha),  levied  each  cutting-cycle  on  harvested  lands, 
changes  the  profit  maximization  objective  function  to: 

An  area  fee  of  $250/ha  is  used  in  the  present  analysis.  An  area  fee  facilitates  government 
rent  capture  relative  to  tree  or  volume  royalties,  though  it  is  not  expected  to  affect 
harvesting  behavior  on  economic  forestlands.  Rather,  those  forest  areas  with  timber 
values  less  than  the  area  fee  are  expected  to  be  taken  out  of  production. 

Per-tree  royalty.  A  per-tree  royalty  of  (r  $/stem)  changes  the  concessionaire's 
objective  function  to 
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The  per-tree  royalty  is  set  at  $70/tree.  This  and  the  other  royalties  are  expected  to  reduce 
the  concessionaire's  profit  margin  on  individual  trees  and  thus  may  lead  to  "high- 
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grading",  or  logging  of  only  those  stems  of  superior  value  and  stem  size  (e.g.,  Repetto 
and  Gillis  1988,  Boscolo  and  Vincent  2000). 

Volume-based  royalty.  Concession  fees  derived  from  a  volume-based  royalty 
(i  $/m3)  would  change  the  objective  function  to 

The  volume-based  royalty  is  set  at  $40/m3  for  the  present  analysis. 

Ad  valorem  royalty.  An  ad  valorem  royalty  is  derived  as  a  percentage  of  the  log 
price  (t  %  log  price)  and  changes  the  logger's  objective  function  to 
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The  ad  valorem  royalty  is  set  at  71%  of  the  log  price  of  harvested  stems. 

Performance  bond  (PB).  Although  uncommon  as  regulatory  instruments  for 
forest  management,  performance  bonds  may  provide  a  useful  incentive  mechanism  for 
logger  compliance  with  desired  management  practices.  For  the  present  study  the 
performance  bond  is  levied  prior  to  management,  and  then  a  proportion  equivalent  to 
logger  compliance  is  refunded  to  the  logger  immediately  after  harvest.  The  performance 
bond  (t  $/ha)  changes  the  objective  function  to 

Compliance  with  the  performance  bond  requires  that  the  residual  merchantable 
stock  population  and  basal  area  of  the  managed  forest  equal  or  exceed  the  residual 
merchantable  stock  of  the  STY  equilibrium  (Eq.  4-13, 4-14).  Any  decrease  in  the 
residual  stock  implies  a  proportional  decrease  in  compliance  with  this  standard  and  thus  a 
proportional  reduction  in  the  performance  bond  refunded  to  the  concessionaire.  The 
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objective  function  is  accordingly  modified  to  reflect  that  profit  maximization  is  defined 
by  the  intensity  of  harvest  of  merchantable  stems  and  the  level  of  compliance  with 
performance  standards  (Eq.  4-12).  In  addition  to  the  forest  growth,  merchantability, 
absolute  harvest,  and  non-negativity  constraints,  the  performance  bond  optimization  is 
constrained  by  criteria  for  compliance  with  sustainability  standards  (Eq.  4-13, 4-14). 
Percent  compliance  is  constrained  to  the  interval  (0,1)  (Eq.  4-15). 
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where 

y  0,  =  Number  of  merchantable  stems  (trees/ha)  in  size  class  I  and  guild./  at  time  t 

f'y  =  Number  of  merchantable  stems  (trees/ha)  in  size  class  i  and  guild./  in  the 
sustained-yield  (SY)  distribution 

p  =  %  compliance  with  SY  restriction 

Bfj-  Mean  basal  area  (mJ)  of  Size  Class  ;'  and  Guild/ 
Intermitting  the  appropriate  fee  level  for  the  performance  bond  is  difficult,  as  setting  the 
bond  too  low  will  have  little  effect  on  management  behavior,  while  setting  it  too  high  will 
discourage  private  investment  in  logging  government-allocated  concessions  (Richards 
2000).  Nonetheless,  for  the  logger  to  have  the  incentive  to  comply  with  management 
standards  and  receive  a  full  refund,  the  bond  must  be  set  sufficiently  high.  Trade  offs 
between  compliance  and  harvest  are  expected  at  lower  performance  bond  levels.  The 
present  PB  optimization  model  allows  for  an  examination  of  "trade-off  solutions",  as  the 
decision  variables  include  both  level  of  compliance  as  well  as  harvest  intensities  per  guild 


and  size  class.  The  performance  bond  level  is  set  at  $750/ha  to  illustrate  one  such  trade- 
off solution  compatible  with  sustained-yield  objectives. 
Results  and  Discussion 

Alternative  regulatory  mechanisms  result  in  changes  in  both  optimal  harvesting 
behavior  and  in  the  distribution  of  timber  rents  between  private  concessionaires  and  the 
government  relative  to  BFL  policies.  The  relative  efficiency  of  regulatory  mechanisms  is 
revealed  by  inspection  of  the  forest  condition,  profitability,  and  allocation  of  stumpage 
between  the  concessionaire  and  the  government.  As  noted  above,  the  relative  efficiency 
of  alternative  policies  is  conditional  on  the  fee  levels  selected  to  meet  the  criterion  of 
maintaining  a  residual  merchantable  timber  stock  approaching  that  of  the  STY  condition. 
Harvesting  Behavior  and  Forest  Condition 

Current  Bolivian  regulatory  policy  (BFL)  comprises  both  technical  management 
standards  (minimum  diameter  limit  and  80%  allowable  cut)  and  an  annual  area  fee  paid 
for  all  forestland  subject  to  logging.  The  annual  area  fee  has  a  negligible  effect  on 
harvest  behavior  for  forestland  with  a  timber  value  exceeding  $5.44/ha,  the  NPV  of 
$l/ha/yr  paid  over  a  40-year  concession  contract  at  a  discount  rate  of  1 8.75%  per  annum. 
At  best,  the  annual  area  fee  provides  an  incentive  to  set  aside  forestland  with  marginal 
timber  value,  which  consists  largely  of  the  inselberg  habitats  with  species  of  no  market 
value  (Guild  3).  The  management  standards  have  the  desirable  effect  of  maintaining 
merchantable  stems  of  smaller  size  classes  (<  45  cm  dbh)  and  maintaining  20%  of  the 
growing  stock  of  merchantable  stems  of  45  cm  dbh  and  greater.  The  residual  distribution 
of  merchantable  stems  under  BFL  deviates  considerably  from  the  STY  condition, 
however,  as  the  optimal  harvest  solution  is  to  extract  all  stems  more  than  50  cm  dbh  and 
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leave  stems  of  the  45  cm  dbh  class  adequate  to  constitute  20%  of  the  growing  stock 

(Table  4-1,  Fig.  4-1). 

The  area  fee  alternative,  which  requires  payment  of  $250/ha  at  each  cutting-cycle 

for  the  area  harvested,  does  not  change  optimal  harvesting  behavior  from  the  liquidation 

harvest  that  would  be  optimal  under  unregulated  management.  The  area  fee  is  paid 

regardless  of  harvest  intensity,  thus  the  concessionaire's  optimal  strategy  for  maximizing 

profits  would  be  to  extract  all  merchantable  timber  at  each  entry.  As  modeled,  the  area 

fee  would  result  in  depletion  of  the  merchantable  stock  over  the  40-year  concession 

period  (Table  4-1,  Fig.  4-1). 

Table  4-1 .  Residual  distribution  of  merchantable  stock  (trees  >  40  cm  dbh  per  ha)  per 
regulatory  policy,  Year  40 


Regulatory 

mechanism 

Dbh  class 

Ad 

(cm) 

STY 

BFL 

Area 

Per-tree 

Volume 

valorem 

PB 

40 

3.2 

6.3 

0.0 

6.3 

3.9 

1.0 

4.6 

45 

2.1 

1.6 

0.0 

3.5 

2.9 

0.8 

3.2 

50 

1.6 

0.0 

0.0 

2.6 

2.0 

0.7 

2.0 

55 

1.3 

0.0 

0.0 

1.3 

1.3 

0.5 

1.1 

60 

1.1 

0.0 

0.0 

0.3 

0.8 

0.3 

0.3 

65 

0.9 

0.0 

0.0 

0.2 

0.5 

0.2 

0.3 

70 

0.7 

0.0 

0.0 

0.0 

0.3 

0.2 

0.2 

75 

0.6 

0.0 

0.0 

0.0 

0.1 

0.1 

0.1 

80 

0.3 

0.0 

0.0 

0.0 

0.1 

0.1 

0.1 

85+ 

0.0 

0.0 

0.0 

0.0 

0.1 

0.1 

0.1 

The  per-tree  royalty  creates  the  incentive  to  leave  smaller  stems  of  lesser  value, 
notably  those  stems  whose  net  value  do  not  exceed  the  $70  per-tree  fee,  such  as  the 
species  of  Guild  1,  which  have  the  lowest  market  value  (Fig.  4-2).  The  per-tree  royalty 
results  in  modest  improvements  toward  achieving  the  STY  condition  relative  to  the  area 
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Figure  4-1 .  Distribution  of  merchantable  population  per  regulatory  policy,  Year  40 
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Figure  4-2.  Distribution  of  merchantable  stems  of  Guild  1  per  regulatory  policy,  Year  40 
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fee  (Table  4-1,  Fig.  4-1),  but  maintains  the  incentive  to  high-grade  the  forest  concession 
and  thus,  as  modeled,  is  not  an  efficient  mechanism  for  attaining  STY  goals. 

The  optimal  harvesting  strategy  under  regulation  by  the  volume-based  royalty 
would  be  to  leave  stems  of  low  value,  notably  those  of  Guilds  1  and  4.  The  merchantable 
stock  of  lesser- valued  species  would  be  maintained  (Fig.  4-1,  4-2),  but  species  whose 
values  exceed  the  volume  fee  ($40/m3)  would  be  depleted  (Fig.  4-3).  The  effect  of  the 
volume  fee  is  similar  to  a  reduction  in  the  price  of  timber.  The  consequence  would  be 
more  selective  high-grading. 
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Figure  4-3.  Distribution  of  merchantable  stems  of  Guild  5  per  regulatory  policy,  Year  40 

The  performance  bond  mechanism  creates  an  incentive  for  partial  to  full 
compliance,  according  to  the  bond  level  (Fig.  4-4).  At  the  $750/ha  bond  level  examined, 
the  concessionaire's  optimal  strategy  would  be  for  approximately  80%  compliance  with 
the  standard  of  maintaining  a  residual  stock  equal  to  that  of  the  STY  condition  (1 1 .7 
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stems/ha)  at  each  cutting-cycle.  In  other  words,  the  optimal  harvest  intensity  is  that 
which  results  in  a  residual  merchantable  population  approximately  20%  lower 
(8.2  stems/ha)  than  that  necessary  for  STY.  Given  the  incentive  for  compliance,  the 
performance  bond  results  in  a  forest  condition  more  closely  approximating  the  STY 
condition  than  the  alternative  regulatory  options  examined  (Table  4-1,  Fig.  4-1).  As 
defined,  however,  the  PB  mechanism  permits  high-grading  of  the  stand  and  satisfaction 
of  performance  criteria  by  leaving  lesser-valued  stems  (Fig.  4-2, 4-3).  In  terms  of 
efficiency,  a  PB  regulatory  policy  is  preferable  to  the  alternative  undifferentiated 
mechanisms,  but  would  be  even  more  effective  with  greater  precision  in  the  definition  of 
performance  standards.  Criteria  for  maintaining  residual  stocking  per  guild  and  size  class 
groupings  would  promote  closer  compliance  with  the  maintenance  of  forest  structure  and 
composition  required  for  STY,  as  defined  by  the  forest  growth  and  optimization  models. 


1.00 


$738.99/ha. 


Performance  bond  (t  $/ha) 

Figure  4-4.  Percent  compliance  (p)  with  the  STY  condition  and  NPV  of  returns  from 
management  ($/ha)  relative  to  performance  bond  level  (t) 
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Marginal  increases  in  compliance  (p)  relative  to  performance  bond  fee  (t)  are  low 
up  to  a  fee  level  of  $650/ha,  at  which  point  compliance  increases  dramatically  with 
increases  in  bond  level  (Figure  4-4).  The  behavior  of  this  compliance  curve  suggests  that 
the  optimal  performance  bond  fee  lies  between  $640/ha  and  $800/ha,  the  precise  level 
defined  by  a  trade-off  between  profitability  and  proximity  to  STY  conditions.  If  the 
marginal  benefit  curve  for  production  of  timber  and  non-timber  benefits  under  STY  were 
known,  the  optimal  fee  level  could  be  more  precisely  estimated  as  that  point  at  which  the 
marginal  cost  of  an  increase  in  fees  equals  to  marginal  benefit  of  greater  compliance. 
The  estimation  of  optimal  fee  levels  for  performance  bond  mechanism  would  be  a 
valuable  extension  of  this  work. 
Profitability  and  Rent  Distribution 

With  the  exception  of  the  area  fee  mechanism,  which  has  no  impact  on  harvesting 
behavior,  regulatory  alternatives  to  the  existing  Bolivian  forestry  law  result  in  lower 
concessionaire  profits  and  greater  government  rent  capture  (Table  4-2). 
Table  4-2.  Profitability  and  rent  distribution  per  regulatory  policy 


Rent 

capture 

(NPV  $/ha) 

Regulatory  policy 

Concessionaire 

Government 

Sustained  timber  yield  (STY) 

257.18 

8.26 

Bolivian  forestry  law  (BFL) 

914.89 

5.33 

Area  fee 

1135.61 

250.26 

Per-tree  royalty 

129.37 

356.32 

Volume-based  royalty 

83.21 

292.13 

Ad  valorem  royalty 

93.98 

1345.17 

Performance  bond  (PB) 

740.13 

148.46 

The  magnitude  of  the  area  fee  has  no  impact  on  harvesting  behavior;  it  simply 
influences  the  net  private  and  governmental  returns  to  management.  Certainly,  if  the  area 
fee  dissipates  all  profits  to  management,  the  concessionaire  is  unlikely  to  log.  As  long  as 
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the  concessionaire  can  expect  profitable  returns  to  management,  however,  the 
concessionaire's  previous  optimal  strategy  of  maximizing  profits  in  liquidating  the 
merchantable  stock  is  unchanged;  the  net  profits  are  just  lower. 

The  royalty  systems  result  in  substantially  lower  private  returns  for  modest 
improvements  toward  the  STY  condition  relative  to  BFL  or,  in  the  ad  valorem  case, 
further  depletion  of  the  forest  stock.  All  alternative  mechanisms  would  result  in  greater 
government  rent  capture  relative  to  current  BFL  conditions. 

The  performance  bond  option  reduces  concessionaire  profits  the  least,  while 
increasing  government  rent  capture  and  resulting  in  greater  proximity  to  the  STY 
condition  relative  to  current  and  alternative  regulatory  mechanisms.  The  challenges  of 
defining  appropriate  performance  criteria  and  enforcing  logger  compliance  under  a 
performance  bond  are  imposing.  Moreover,  high  performance  bond  fees  are  required  to 
induce  favorable  levels  of  compliance  with  desired  forest  management  standards. 
Excessive  fees  will  discourage  private  investment  in  logging  government-allocated 
concessions.  The  incentive  for  compliance  created  by  a  performance  bond  mechanism  is 
particularly  compelling,  though  the  complexity  of  its  application  may  be  limiting. 

Conclusions 

Alternatives  to  current  Bolivian  regulatory  policy  for  forest  concessions  may 
provide  more  effective  fiscal  incentives  for  sustainable  forestry,  while  improving 
government  rent  capture.  The  goal  of  changing  logger  behavior  toward  greater 
compliance  with  STY  standards  may  be  achieved  by  alternative  royalty  systems  or  a 
performance  bond  mechanism,  given  appropriate  fee  levels.  An  area  fee,  unless 
complemented  by  technical  standards  for  management  such  as  the  Normas  Tecnicas,  is 
unlikely  to  change  the  concessionaire's  incentive  to  liquidate  all  merchantable  timber  in 
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initial  harvest  entries.  For  all  alternative  regulatory  mechanism  forms  examined  in  this 
study,  optimal  concessionaire  behavior  is  to  high-grade  the  forest,  as  revealed  in 
examination  of  the  residual  stock  of  stems  of  the  highest  value,  which  constitute  Guild  5 
(Fig.  4-3).  This  optimal  "high-grading"  behavior  may  be  modified  with  differentiation  of 
regulatory  fees  (Amacher  et  al.  2001),  such  as  higher  per-tree  and  volume  royalties  for 
more  valuable  relative  to  lesser-valued  species  and  more  precise  STY  compliance 
standards  for  the  PB  mechanism.  Presently,  the  policy  alternatives  are  defined  as 
undifferentiated  mechanisms,  meaning  that  all  stems  are  subjected  to  the  same  level  of 
fiscal  regulation.  This  regulatory  approach  is  considerably  more  straightforward  to 
implement,  but  does  not  result  in  the  differential  harvesting  behavior  that  would  be 
necessary  to  compel  loggers  to  conserve  the  growing  stock  of  more  valuable  species. 
Examining  better  differentiated  forms  of  the  performance  bond  and  royalty  mechanisms 
would  be  a  valuable  extension  of  this  research. 

Establishing  straightforward  and  meaningful  standards  by  which  the  appropriate  fee 
levels  may  be  estimated  and  concessionaire  management  performance  may  be  evaluated 
is  problematic  because  of  the  complexity  of  mixed  tropical  forests  and  the  difficulty  of 
estimating  forest  conditions  necessary  for  sustained  timber  productivity.  Appropriate 
regulatory  mechanisms  must  balance  efficiency  gains  of  improved  compliance  with  STY 
standards  and  desirable  distribution  of  resource  rents  with  the  practical  aspects  of 
implementation.  Regulatory  mechanisms  that  are  too  complex  will  likely  be  poorly 
implemented  and  are  unlikely  to  achieve  the  theoretical  efficiency  gains  for  which  they 
are  designed. 
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Moreover,  the  design  of  fiscal  regulatory  policies  has  implications  across  a  broader 
landscape  of  land  uses.  While  higher  forest  fees  seem  necessary  to  counter  incentives  for 
high-grading  and  to  constrain  optimal  behavior  to  the  conditions  necessary  for 
sustainability,  increasing  forest  fees  lowers  the  profits  of  forest  management  for  timber. 
A  decrease  in  the  profitability  of  logging  implies  that  natural  forests  of  marginal 
economic  value  will  drop  out  of  timber  production.  Higher  taxation  reduces  the  incentive 
for  maintaining  forestland  in  timber  production  for  all  but  the  most  highly  valued  species 
and  creates  greater  incentive  to  convert  natural  forests  to  plantations  or  non-forest  uses. 
The  optimal  regulatory  policy  is  likely  a  compromise  between  efficient  fiscal 
mechanisms  and  modest  fee  increases,  complemented  by  the  enforcement  of  appropriate 
technical  standards,  or  best  management  practices. 


CHAPTER  5 
CONCLUSIONS 

The  present  study  of  optimal  management  regimes  for  a  Chiquitano  tropical  dry 
forest  offers  an  indication  of  the  biological  and  economic  requirements  for  sustained 
timber  yield  (STY)  and  sustainable  forest  management  (SFM),  given  a  land  use  objective 
of  profit  maximization  in  logging  and  the  biological,  market  and  regulatory  constraints  on 
its  achievement.  Results  generated  by  the  forest  growth  model  and  optimization  analyses 
are  conditional  on  the  modeling  assumptions  required  for  a  tractable  analysis.  The 
accuracy  of  these  bioeconomic  models  fundamentally  determines  the  validity  of  estimates 
concerning  the  conditions  required  for  sustainable  management  and  simulations  of 
management  impacts  on  the  natural  forests.  Constraining  the  harvest  horizon  to  a  40-year 
concession  period  lends  greater  integrity  to  the  analysis  of  logging  behavior,  but  does  not 
resolve  all  issues  caused  by  data  limitations  and  restrictive  modeling  assumptions. 

These  limitations  notwithstanding,  the  study  offers  several  useful  contributions  to 
research  on  the  biological  and  economic  requirements  for  sustainable  forest  management. 
The  novel  application  of  the  multinomial  logit  (MNL)  regression  method  for  estimating 
matrix  transition  models  demonstrates  important  benefits  relative  to  previous  techniques 
employed  in  matrix  model  estimation.  The  MNL  method  corrects  for  variance  in  the  data 
and  model  estimation  errors  imposed  by  limited  samples.  Moreover,  the  MNL  form 
offers  compatible  simulation  alternatives,  including  deterministic  and  stochastic 
projection  in  both  static  and  dynamic  model  forms.  The  study  demonstrates  that 
sustaining  the  productive  capacity  of  the  Chiquitano  forest  for  timber  can  be  compatible 
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with  generating  positive  financial  returns  from  management  and  with  improved 
ecosystem  conservation.  The  integrated  bioeconomic  optimization  model  reveals  the 
governing  influence  of  forest  growth,  cutting-cycle  length,  and  timber  merchantability  on 
the  sustainable  management  solutions.  Finally,  the  analysis  of  regulatory  policies 
demonstrates  the  effects  and  relative  efficiencies  of  alternative  fiscal  regulatory 
instruments  as  means  of  inducing  logger  compliance  with  desirable  forest  management 
practices. 

The  extraordinary  discount  rate  of  1 8.75%  applied  in  this  study  and  the  slow 
growth  rates  of  merchantable  timber  species  suggest  that  short  cutting-cycles  and 
measured,  selective  harvesting  promise  the  highest  financial  returns  to  sustainable 
management  of  the  Chiquitano  forest.  This  high  discount  rate  also  substantially  inflates 
the  opportunity  costs  of  more  sustainable  management  relative  to  unconstrained 
harvesting.  Given  the  low  value  growth  rate  of  tropical  hardwoods,  however,  the 
requirement  for  leaving  merchantable  timber  in  the  forest  imposed  by  both  the 
sustainable  forestry  scenarios  (STY  and  SFM)  implies  that  these  options  will  consistently 
generate  inferior  returns  to  unconstrained  liquidation  of  merchantable  timber  in  initial 
harvest  entries. 

Results  of  the  SFM  optimization  indicate  that  maintaining  forest  structure  and 
composition  can  be  compatible  with  producing  sustained  timber  yield  to  meet  the 
economic  objective  of  generating  profitable  returns  from  management.  Moreover,  for  the 
particular  case  of  Chiquitano  tropical  dry  forests,  both  STY  and  SFM  goals  appear 
compatible  with  the  maintenance  of  forest  wildlife  values.  More  information  on  the 
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values  of  timber  species  for  wildlife  taxa  other  than  vertebrate  frugivores  would  allow  for 
a  more  conclusive  assertion  of  this  apparent  compatibility. 

The  finding  that  the  opportunity  costs  of  sustainable  management  alternatives  are 
high,  relative  to  unsustainable  logging  methods,  should  not  imply  that  the  pursuit  of 
sustainability  is  an  unworthy,  or  impractical,  objective.  The  private  opportunity  costs  of 
sustainable  management  may  be  high,  but  the  public  opportunity  costs  of  the  progressive 
depletion  of  natural  forest  timber  stocks  and  concomitant  degradation  of  non-timber 
benefits  produced  by  natural  forests  that  are  expected  to  occur  under  conventional, 
unsustainable  logging  are  arguably  much  higher.  The  critical  issues  for  this  and  similar 
research  on  sustainable  tropical  forestry  are  whether  sustainable  management  for  timber 
and  non-timber  benefits  is  biologically  feasible  and  economically  viable.  If  so,  the  goal 
of  establishing  permanent  production  forests  to  meet  demands  for  tropical  hardwoods  is 
both  a  worthy  and  practical  alternative  to  extensive  high-grading  and  degradation  of 
natural  forests  for  short-term  timber  production. 

Examination  of  the  economic  and  biological  impacts  of  optimal  management  under 
constraints  imposed  by  current  Bolivian  forestry  regulations  (BFL)  indicates  that, 
although  the  financial  returns  to  management  are  appreciable,  BFL  standards  are  not 
likely  to  achieve  sustainability  goals.  This  conclusion  cannot  be  drawn  without  caveats, 
given  the  simplifying  assumptions  accepted  in  analyzing  the  BFL  scenario  with  the 
available  forest  growth  and  optimization  models.  Optimization  results  indicate  that  tree 
diameter  and  allowable  cut  restrictions  of  BFL  provide  the  logger  with  an  incentive  to 
deplete  large  stems  of  high  value,  a  practice  unlikely  to  result  in  the  maintenance  of  forest 
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structure  and  composition  necessary  to  achieve  sustained  timber  yield  by  present  model 
estimates. 

Alternatives  to  current  Bolivian  regulatory  policy  for  forest  concessions  may 
provide  more  effective  fiscal  incentives  for  sustainable  forestry,  while  improving 
government  rent  capture.  The  goal  of  compelling  logger  behavior  toward  greater 
compliance  with  sustainability  standards  may  be  more  efficiently  achieved  by  alternative 
royalty  systems,  or  a  performance  bond  mechanism,  given  appropriate  fee  levels.  High- 
grading  is  expected  to  result  from  greater  fiscal  regulation,  however,  if  such  regulations 
are  undifferentiated  across  timber  value  classes.  Examining  more  precise,  differentiated 
forms  of  performance  bond,  and  royalty  mechanisms,  would  be  a  valuable  extension  of 
this  research.  Still,  establishing  meaningful  standards  by  which  management 
performance  may  be  evaluated  is  problematic  because  of  the  complexity  of  mixed 
tropical  forests  and  the  difficulty  of  estimating  forest  conditions  necessary  for  sustained 
timber  productivity.  Appropriate  regulatory  mechanisms  must  balance  efficiency  gains 
of  improved  compliance  with  sustainability  standards  and  of  a  more  desirable  distribution 
of  resource  rents  with  the  practical  aspects  of  their  implementation  and  enforcement. 

Numerous  interesting  extensions  of  this  work  are  possible.  The  assumption  of 
stationarity  for  transition  probabilities  of  the  matrix  model  limits  the  examination  of 
forest  dynamics.  An  inspection  of  forest  impacts  and  sustainability  conditions  under 
dynamic  modeling  would  make  better  use  of  the  research  potential  of  the  MNL  model. 
Recalculation  of  transition  probabilities  per  harvest  episode  and  silvicultural  treatment 
may  offer  greater  insight  into  the  expected  forest  response  to  management  and  parameters 
for  optimal  harvest.  The  optimal  path  of  regulated  harvest  necessary  to  maximize  profits 
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while  converting  the  forest  to  a  sustained-yield  state  is  another  compelling  extension. 
The  period  of  regulation  required  to  establish  STY  conditions  and  the  financial 
implications  of  this  regulation  may  offer  important  insights  into  the  management  needs 
for  establishing  a  regulated  forest.  The  integrated  optimization  model  enables  an 
examination  of  conservation  and  management  alternatives  proposed  in  recent  conceptual 
studies.  Examples  include  the  "log  and  protect  model"  suggested  by  Rice  et  al.  (2001), 
the  effects  of  a  price  premium  for  certified  timber  with  and  without  effective  harvest 
regulation,  and  the  effects  of  a  market  for  carbon  sequestration. 

Results  of  the  present  research  offer  a  modest,  but  hopefully  useful,  contribution  to 
scientific  understanding  of  the  complex  and  contentious  issue  of  sustainable  forestry. 
Improving  the  precision  and  dynamics  of  the  modeling  tools  developed  in  this  study  will 
enable  further  and  more  practical  contributions  to  management,  economic  analysis,  and 
policy  design. 


APPENDIX  A 
GUILD  CLASSIFICATION  OF  LOMERIO  TREE  SPECIES 


Forest 

Regen.  shade  Mature  shade    Timber 

Wildlife 

Scientific  name            Guild  habitat* 

tol.  (1-3)       tol.  (0,1)     value  (0-3) 

value  (1-3) 

Allophylus  edulis 

6 

3                   1                 0 

1 

Ampelocera  ruizii 

I          7 

3                   1                 0 

2 

Aspidosperma  cylindrocarpon 

I           7 

2                   1                  1 

1 

Aspidosperma  pyrifolium 

I          7 

2                   1                  1 

1 

Attalea  phalerata 

I           8 

2                   1                 0 

3 

Bougainvillea  sp. 

I          5 

2                   1                 0 

1 

Bouganvillea  modesta 

I          5 

2                   1                 0 

1 

Calyptranthes  sp. 

1          6 

3                   1                 0 

2 

Campomanesia  aromatica 

I          7 

2                   1                 0 

3 

Capparis  prisca 

I          7 

3                   1                 0 

3 

Cariniana  ianeirensis 

I          7 

2                   1                  1 

2 

Chrysophyllun  gonocarpum 

I          7 

3                   1                 0 

3 

Dalbergia  riparia 

1          8 

2                   1                 0 

1 

Duguettia  guitarensis 

I          8 

3                   1                 0 

3 

Erythroxylum  sp. 

I          6 

3                   1                 0 

1 

Esenbeckia  almcmillea 

I          6 

3                   1                 0 

1 

Galipea  trifoliata 

I          6 

3                   1                 0 

1 

Gallesia  integrifolia 

I          5 

3                   1                  1 

1 

Genipa  americana 

I          7 

2                  1                 0 

3 

Inga  marginata 

1          7 

2                   1                 0 

3 

Lonchocarpus  guillemineanus 

1          7 

2                  0                 1 

1 

Myrciaria  cauliflora 

I          8 

3                   1                 0 

3 

Myrtaceae  sp. 

[          7 

2                   1                 0 

3 

Neea  hermaphrodita 

[          6 

3                   1                 0 

2 

Phyllostylon  rhamnoides 

I          7 

3                   1                 0 

1 

Piper  alboreum 

8 

3                   1                 0 

2 

Qualea  acuminata 

1          7 

2                   1                 0 

1 

Rhamnidium  elaeocarpum 

I          7 

2                   1                 0 

3 

Rheedia  acuminata 

8 

3                  2                 0 

3 

Salacia  elliptica 

8 

3                   1                 0 

3 

Sorocea  saxicola 

7 

3                   1                  1 

3 

Syagrus  sancona 

8 

3                   1                 0 

3 
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Forest  Regen.shade  Mature  shade    Timber      Wildlife 
Scientific  name  Guild  habitat*    tol.  (1-3)       tol.  (0,1)    value  (0-3)  value  (1-3) 


Talisia  esculenta 

7 

3 

1 

0 

3 

Trichilia  palliata 

6 

3 

2 

0 

1 

Triplaris  americana 

8 

2 

1 

0 

1 

Vitex  cymosa 

7 

3 

1 

0 

3 

Ximenia  americana 

8 

2 

1 

0 

3 

Cariniana  domestica 

2 

7 

0 

2 

2 

Cecropia  concolor 

2 

8 

0 

0 

3 

Cyclolobium  blanchetianum 

2 

7 

0 

0 

1 

Ficus  gomelleira 

2 

7 

0 

0 

3 

Guazuma  ulmifolia 

2 

8 

1 

1 

3 

Hymenea  courbaril 

2 

7 

0 

2 

3 

Samanea  tubulosa 

2 

8 

1 

0 

3 

Sapindus  saponaria 

2 

8 

1 

0 

3 

Sapium  longifolium 

2 

8 

1 

0 

1 

Tabebuia  serratifolia 

2 

7 

0 

2 

2 

Zeyheria  turberculosa 

2 

7 

0 

1 

2 

Agonanda  brasilensis 

3 

3 

1 

0 

2 

Bauhinia  ungulata 

3 

2 

0 

0 

Celtis  iguanea 

3 

2 

1 

0 

Celtis  spinosa 

3 

2 

0 

0 

Cereus  spp. 

3 

3 

2 

0 

0 

Cochlospermum  vitifolium 

3 

2 

0 

0 

Eriotheca  roseorum 

3 

3 

0 

0 

Heliocarpus  americanus 

3 

1 

1 

0 

Jacaranda  cuspidifolia 

3 

3 

0 

0 

Luehea  peniculata 

3 

3 

1 

0 

Pseudobombax  marginatum 

3 

3 

0 

0 

Schefflera  morototoni 

3 

2 

0 

0 

Sebastiana  brasiliensis 

3 

1 

1 

0 

Stryphnodendron  guianense 

3 

3 

0 

0 

Tabernaemontana  sp. 

3 

4 

1 

0 

Trigonia  boliviano 

3 

2 

0 

0 

Urera  baccifera 

3 

1 

0 

0 

Vocnysia  mapirensis 

3 

2 

0 

0 

Amburana  cearensis 

4 

5 

0 

3 

Anadenanthera  colubrina 

4 

5 

0 

2 

Astronium  urundeuva 

4 

5 

0 

2 

Cedrelafissilis 

4 

5 

0 

3 

Ceiba  samauma 

4 

5 

0 

0 

2 
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Forest 

Regen.shade Mature  shade 

Timber 

Wildlife 

Scientific  name 

Guild  habitat* 

tol.  (1-3) 

tol.  (0,1) 

value  (0-3) 

value  (1-3) 

Centrolobium  microchaete 

4 

4 

1 

0 

2 

1 

Chorisia  speciosa 

4 

5 

1 

0 

0 

2 

Cordia  glabrata 

4 

5 

1 

0 

3 

Mimosa  sp. 

4 

4 

1 

0 

0 

Physocalymma  scaberhmum 

4 

5 

2 

0 

0 

Piptadenia  viridifolia 

4 

5 

2 

0 

0 

Platymiscium  ulei 

4 

5 

1 

0 

0 

Pterogyne  nitens 

4 

5 

1 

0 

2 

Schinopsis  brasilensis 

4 

5 

1 

0 

3 

Spondias  mombin 

4 

5 

2 

0 

0 

Terminalia  oblonga 

4 

5 

1 

0 

1 

Acacia  albocorticata 

5 

5 

1 

0 

0 

2 

Acosmiun  cardenasii 

5 

5 

3 

1 

0 

Aspidosperma  nobile 

5 

5 

2 

0 

1 

Aspidosperma  rigidum 

5 

5 

2 

0 

2 

Caesalpinia  paraguarensis 

5 

4 

1 

0 

0 

Caesalpinia  pluviosa 

5 

5 

2 

0 

0 

2 

Capparis  retusa 

5 

4 

3 

1 

0 

3 

Casearia  gossipiosperma 

5 

4 

3 

1 

0 

3 

Combretum  leprosum 

5 

5 

2 

1 

0 

1 

Copaifera  sp. 

5 

5 

3 

1 

0 

2 

Dictyoloma  peruviana 

5 

5 

1 

0 

0 

1 

Enterolobium  contortisiliquum 

5 

5 

1 

0 

0 

3 

Guibourtia  chodatiana 

5 

5 

3 

1 

1 

3 

Jatropha  minuscula 

5 

5 

1 

0 

0 

1 

Machaerium  aculifolium 

5 

5 

3 

1 

0 

1 

Machaerium  hirtum 

5 

5 

1 

0 

0 

1 

Machaerium  jacarandifolium 

5 

4 

3 

0 

0 

1 

Machaerium  scleroxylon 

5 

5 

3 

0 

3 

1 

Madura  tinctoria 

5 

5 

1 

0 

0 

2 

Opuntia  brasiliensis 

5 

5 

3 

0 

0 

3 

Pereskia  sacharosa 

5 

5 

2 

0 

0 

1 

Platypodium  elegans 

5 

5 

3 

0 

1 

1 

Pogonopus  tubulosus 

5 

5 

3 

1 

0 

1 

Rollinia  emarginata 

5 

4 

3 

1 

0 

3 

Simira  rubescens 

5 

5 

3 

0 

0 

1 

Swartzia  jorori 

5 

4 

3 

1 

0 

3 

Zanthoxylum  sp. 

5 

5 

1 

0 

0 

2 

K4 


*  Lomerio  forest  habitats  (T.  Fredericksen.  personal  comment) 

1 :  Short-lived  (3-4  years),  early  successional 

2:  Long-lived  colonizers  of  disturbed  sites 

3:  Borders  of  rock  outcrops 

4:  Generalists 

5:  Mature  upland  forest  canopy 

6:  Mature  upland  forest  understory 

7:  Riparian  forest  canopy 

8:  Riparian  forest  understory 


85 


APPENDIX  B 
FOREST  HARVESTING  COSTS  PER  COMPONENT 


Direct  costs 

$US/ha 

$US/m3 

Road  construction  and  maintenance 

18.81 

3.33 

Management  plan  (compounded  2  yrs  @  18.75%) 

0.91 

0.26 

Inventory 

0.31 

0.06 

Plan  (data  analysis,  mapping,  write-up) 

0.36 

0.09 

Census  of  merchantable  stems 

10.04 

2.49 

(compounded  1  year  @1 8.75%) 

Annual  Operating  Plan 

0.86 

0.10 

(compounded  1  year  @18.75%) 

Felling  &  bucking 

Labor  (sawyer  +  assistant) 

2.57 

0.49 

Chainsaw(Stihl070) 

Depreciation 

0.24 

0.05 

Chains 

0.20 

0.04 

Blades 

0.14 

0.03 

Pinions 

0.03 

0.01 

Files 

0.03 

0.00 

Spark  plugs 

0.01 

0.00 

Filters 

0.02 

0.00 

Fuel 

0.16 

0.02 

Motor  oil 

0.06 

0.01 

Safety  equipment 

0.01 

0.00 

Skidding 

Labor  (skidder  operator  +  assistant) 

3.08 

0.59 

Skidder(CAT518) 

Depreciation 

4.74 

0.96 

Cables 

0.33 

0.06 

Tires 

0.34 

0.08 

Filters 

0.03 

0.00 

Annual  maintenance 

1.01 

0.19 

Fuel 

2.61 

0.48 

Motor  oil 

0.35 

0.05 

Hydrolic  fluid 

0.06 

0.01 

Safety  equipment 

0.01 

0.00 
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Measuring  and  bucking  at  log  landing 

Labor  (sawyer  +  technician) 

6.21 

1.18 

Radio 

0.08 

0.01 

Fuel 

2.61 

0.48 

Motor  oil 

0.35 

0.05 

Loading 

Labor  (loader  +  assistant) 

2.91 

0.55 

Loader  (CAT  950B) 

Depreciation 

8.18 

1.12 

Tires 

0.54 

0.07 

Filters 

0.03 

0.00 

Annual  maintenance 

1.26 

0.17 

Fuel 

3.89 

0.53 

Motor  oil 

0.23 

0.03 

Hydrolic  fluid 

0.49 

0.07 

Supervision  of  harvest  operations 

9.76 

1.85 

(activities  of  forest  engineer  and  site  manager) 

Transport  (forest  to  Conception  mill,  -30  km.) 

33.70 

5.72 

Indirect  costs 

$US/ha 

$US/m3 

Taxes  and  fees 

Area  fee  (Patente) 

1.00 

6.35 

Registration:  Camara  Forestal  (CF) 

0.10 

0.02 

Registration:  Superintendencia  Forestal 

0.29 

0.04 

Registration:  Export  Chamber 

0.15 

0.02 

Forest  Origin  Certificate  (CFO) 

1.46 

0.28 

Aggregate  Valor  Tax  (IV A) 

0.27 

0.04 

Training 

1.24 

0.20 

Overhead  (10%  of  direct  costs) 

11.69 

2.11 

Employee  benefits  (included  in  labor  cost  calculation) 

Pension  fimd  (AFP) 

2.5% 

Medical  insurance 

10.0% 

Vacation  (days/yr) 

17.50 

Bonus  (extra  month  pay) 

1.00 

Total  costs 

$US/ha 

$US/m3 

Variable  Costs 

12.11 

Fixed  Costs 

62.51 

Annual  Costs 

1.55 
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APPENDIX  C 
MARKET  VALUE  OF  MERCHANTABLE  TIMBER  SPECIES 


Mean  sawnwood  price 

Mean  roundwood  pricea'b 

(FOB  Santa  Cruz') 

(FOB  Forest  Mill} 

Scientific  name 

B^f 

SUSftf 

Bs/m3 

$US/m3 

Amburana  cearensis 

3.50 

0.52 

601.44 

89.90 

Anadenanthera  colubrina 

2.40 

0.36 

344.92 

51.56 

Aspidosperma  nobile 
Aspidosperma  spp. 
Astronium  urundeuva 

2.75 
2.65 
2.65 

0.41 
0.40 
0.40 

426.54 
403.22 
403.22 

63.76 
60.27 
60.27 

Cariniana  domestica 

2.00 

0.30 

251.64 

37.61 

Cariniana  ianeirensis 

1.75 

0.26 

193.34 

28.90 

Cedrelafissilis 
Centrolobium  microchaete 

4.00 
3.00 

0.60 
0.45 

718.04 
484.84 

107.33 
72.47 

Cordia  glabrata,  C.  alliodora 

3.80 

0.57 

671.40 

100.36 

Gallesia  integrifolia 
Guibourtia  chodatiana 

1.55 
3.00 

0.23 
0.45 

146.70 
484.84 

21.93 
72.47 

Hymenea  courbaril 

2.90 

0.43 

461.52 

68.99 

Lonchocarpus  guillemineanus 
Machaerium  scleroxylon 
Phyllostylon  rhamnoides 

1.30 

10.51 

2.40 

0.19 

1.57 
0.36 

88.40 

2237.01 

344.92 

13.21 

334.38 

51.56 

Schinopsis  brasilensis 
Sorocea  saxicola 

3.75 
2.40 

0.56 
0.36 

659.74 
344.92 

98.62 
51.56 

Sweetiafruticosa 

2.90 

0.43 

461.52 

68.99 

Tabebuia  serratifolia 

2.75 

0.41 

426.54 

63.76 

Terminalia  oblonga 

2.65 

0.40 

403.22 

60.27 

a.  Less  transport  (0.36  Bsfof)  and  milling 

b.  233.2  bf/m3  roundwood 

costs  (19.38  Bs/m3) 

APPENDIX  D 
DEFINITION  OF  OPTIMIZATION  VARIABLES 

yt  [>"!><]>  number  of  trees  in  size  class  ;',  guild/  at  time  / 

y"  \yi/],  undisturbed  steady-state  or  "climax"  equilibrium  distribution  of  trees 

y<  [  J!rt™]>  number  of  trees  in  size  class  i,  guild/  at  time  I 

y  [y/]>  number  of  trees  in  size  class  i,  guild/  in  the  STY  steady-state 

y"  [y/j"],  number  of  trees  in  size  ;',  guild/  in  the  SFM  steady-state 

h,  [hiji],  number  of  stems  harvested  in  size  class  /',  guild/  at  time  t 

h*  [htj'],  number  of  stems  harvested  in  size  class  z,  guild/  in  the  STY  steady-state 

h"  [htj"},  number  of  stems  harvested  in  size  class  ;',  guild/  in  the  SFM  steady-state 

G,         Matrix  of  transition  probabilities  for  trees  of  guild/,  which  define  the  movement 
of  trees  into  size  classes  (i)  during  a  5-year  growth  period  (6) 

r,  [r/J,  number  of  recruits  into  the  10  cm  size  class  of  guild/  during  the  5-year 

growth  interval  r  to  t+8 

Z         Objective  function  of  the  optimization 

t  =  y8   Cutting-cycle  length,  defined  as  the  number  (y)  of  5-year  growth  periods  (0) 

T       Harvest  horizon 

d       Discount  rate 

Bij      Mean  basal  area  (m  )  of  a  stem  in  size  class  i  and  guild/ 

cq,       Merchantable  proportion  of  stems  in  guild/ 

Sj       Proportion  of  merchantable  species  in  guild/ 

fj         Proportion  of  stems  of  merchantable  form  in  guild/ 

Pj         Market  value  of  merchantable  species  ($/m3)  in  guild/ 
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vij         (Pj  "  Qw> net  value  of  an  extracted  stem  of  size  i  in  guild/  in  200 1  $US/m 

C         Variable  costs  ($/m  )  of  felling,  skidding  and  log  deck  operations 

F         Fixed  costs  ($/ha)  of  planning  and  capital  costs  including  depreciation 

A         Annual  costs  ($/ha)  including  an  area  fee  of  $  1  /ha  for  the  concession  and 
inscription  fees  to  government  agencies  and  trade  associations 

iff,       Mean  extracted  volume  of  a  tree,  estimated  as  85%  of  mean  volume  per  size  class 
(0  and  guild  (/),  assuming  15%  of  merchantable  volume  lost  in  logging  waste 

Vj,        Net  value  ($/m3)  of  an  extracted  stem  of  size  i  in  guild/ 

5  Annual  concession  fee  or  "patente"  ($l/ha)  per  Bolivian  forestry  law 

t  Regulatory  fee  defines  per  fiscal  instrument  (see  Appendix  F) 

p  %  compliance  with  STY  restriction  measured  for  the  performance  bond 

d         [dy],  vector  of  negative  deviations  of  the  number  of  trees  of  size  i  and  guild/  in 
the  yt  from  the  climax  state  y\  expressed  as  an  absolute  value 

a         [di/],  vector  of  positive  deviations  of  the  number  of  trees  of  size  i  and  guild/  in  yt 
from  the  climax  state  y 
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APPENDIX  E 
MATHEMATICAL  DESCRIPTION  OF  MANAGEMENT  ALTERNATIVES 

Unconstrained  Harvest  (U) 


r-l 

I  .0 


h,  +  y,  -  C  ( y,  _,  -  h, ., )  =  J]  G V  Forest  growth 


*(i  ~wjy,jt  -  0  Merchantable  stock,  initial  cutting- 

cycle  (CC) 
r-i 
fy»  -  C/V-i  ~  X  mjG%-}  £  °  Merchantable  stock,  all  ensuing  CC 

4=0 
35 

2  fy ■  =  0  Merchantable  dbh  limit 

i-O 

yi-h,>0  Absolute  harvest  level 

h,  >0,y,>0  Non-negativity 


Bolivian  Forestry  Law  (BFL) 

&v'h,-F     /i_(i  +  <y)-n 


h,+y,  -Gr(y,_t  -h,_,)  =  ^Gkr  Forest  growth 


h,,t  ~  co^ij,  <  0  Merchantable  stock,  initial  CC 

M 

^»  -  GYii-i  -  £  cop''r„,-\  s  0  Merchantable  stock,  all  ensuing  CC 

*-o 


91 


X  E  Z  K  =  °  Lesal  dbh  limit 

,.o    ;     i 

85  85 

£  Ap  -0.8a>y  Y/"*  S  0  Allowable  cut 

MS  (=45 

y,-h,  >  0  Absolute  harvest  level 

h,  >0,y,>0  Non-negativity 


Sustained  Timber  Yield  (STY) 


■^Vh.-F     J\-(\  +  syT) 

maxZlh  „,  =  V — ! --A\  — * '—  ,  s.t. 


Ch'  +  [I-Gr  )y  =  £  Gkr  Sustained-yield  forest  growth 

fr-0 

K   ~a,,y,i   -0  Merchantable  stock 

35 

vVt   =0  Merchantable  dbh  limit 

M    j 

y*-h*>0  Absolute  harvest  level 

h*  >  0,  y*  >0  Non-negativity 

Sustainable  Forest  Management  (SFM) 

i    / 

vft    -FJ_>0  Non-negative  NPV 

(l  +  <J)'fl      <5 

,  .         r-i 

Gr/T  +  (/  -  G'  )y"  =  Y,  G V  Sustained-yield  forest  growth 
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hj    ~0>j)'ij    ^0  Merchantable  stock 

35 

2  X  V  =  "  Merchantable  dbh  limit 

y**  -<t  +  d  =y  Integrity  of  deviation  calculation 

y**  -  h**  >  0  Absolute  harvest  level 

h**,  y**,  <t,  d  >  0  Non-negativity 
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APPENDIX  F 

MATHEMATICAL  DESCRIPTION  OF  OPTIMIZATION  SCENARIOS  WITH 

REGULATORY  MECHANISMS 

Sustained  Timber  Yield  Management  (STY) 


r^-C^-Fh-^Sr) 


max  Zlh  „,  =  7   f  y  - 

i      ]    r-0  (l  +  o)  ^  O 

G'ti  +  \I-G'  )y  =Y_,Gkr  Sustained-yield  forest  growth 

htJ   -e>,yt   ^0  Merchantable  stock 

35 

/, 7,fta    =0  Merchantable  dbh  limit 

1-0     j 

y*-h*>0  Absolute  harvest  level 

h*  >  0,  y*  >0  Non-negativity 

Bolivian  Forestry  Law  (BFL) 

where,  c,  =  annual  area  fee  "patente"  of  ($l/ha),  s.t. 

r-i 
h,+y,-Gr(yl_,-hl_l)  =  YJGkr  Forestgrowth 

i-0 

hljt  -  (Ojytj,  <  0  Merchantable  stock,  initial  cutting 

cycle  (CC) 
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h¥  -  Grym,j,-\  -  ^  ajGkrjjl_l  <  0  Merchantable  stock,  all  ensuing  CC 

40 

S  X  Z  K  -  °  Le§al  dbh  limit 

1-0     j       , 

85  85 

£  \,  -  0.8oy  £  v"/,,  <  0  Allowable  cut 

(=45  (=45 

yt-h,  >  0  Absolute  harvest  level 

hi  >0,y,>0  Non-negativity 


Area  Fee  (t  $/ha) 


r-l 

/j;  +;',  -  C  (x_,  -/!,_,)  =  ^  G*r  Forest  growth 

i-0 

\i  ~a/y^  -  0  Merchantable  stock,  initial  CC 

r-i 
Jfy,  -  G'yVi  -  ^  a>JGkrl]l_l  <  0  Merchantable  stock,  all  ensuing  CC 

i-0 

35 

2j  fy  =  0  Merchantable  dbh  limit 

1-0 

yt-h,>0  Absolute  harvest  level 

hi  >0,y,>0  Non-negativity 


Note:    All  subsequent  regulatory  mechanisms  are  subject  to  the  constraints  defined 
for  the  area  fee  mechanism. 
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Per-Tree  Royalty  (t  $/stem) 


Volume-Based  Royalty  (t  $/m  ) 


y4rf  valorem  Royalty  (t  %  log  price): 

(M)  rrs     (\+sy         {     s     ) 

Performance  Bond  (x  $/ha) 

constraints  for  "Area  Fee"  (above)  and 
8s       r 


1-40    )     (-0 
(=40    7     i-0 

P<1 

/3>0 

where,    J^"  number  of  merchantable  stems  (trees/ha)  in  size  class  i 
and  guild/  at  time  t 

jP  y  "  number  of  merchantable  stems  (trees/ha)  in  size  class  i 
and  guild/  in  the  STY  distribution 

p   =    %  compliance  with  STY  restriction 

Bjj=     mean  basal  area  (m2)of  a  stem  in  size  class  i  and  guild/ 
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